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Spectral hardening as a viable alternative to disc 
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ABSTRACT 

Constraining the accretion flow geometry of black hole binaries in outburst is com- 
plicated by the inability of simplified multi-colour disc models to distinguish between 
changes in the inner disc radius and alterations to the emergent spectrum, parame- 
terised by the phcnomcnological colour correction factor, / co i. We analyse Rossi X-ray 
Timing Explorer observations of the low mass Galactic black hole X-ray binary, GX 
339-4, taken over seven epochs when the source was experiencing a state transition. 
The accretion disc component is isolated using a pipeline resulting in robust detections 
for disc luminosities, 10~ 3 < -^disc/^Edd 0-5. Assuming that the inner disc remains 
situated at the innermost stable circular orbit over the course of a state transition, we 
measure the relative degree of change in f co \ required to explain the spectral evolution 
of the disc component. A variable / co i that increases by a factor of ~ 2.0 — 3.5 as 
the source transitions from the high/soft state to the low/hard state can adequately 
explain the observed disc spectral evolution. For the observations dominated by a 
disc component, the familiar scaling between the disc luminosity and effective tem- 
perature, Ldisc x ^efff' i s observed; however, significant deviations from this relation 
appear when GX 339-4 is in the hard intermediate and low/hard states. Allowing for 
an evolving f co \ between spectral states, the £disc — law is recovered over the full 
range of disc luminosities, although this depends heavily on the physically conceivable 
range of f co \. We demonstrate that physically reasonable changes in / co i provide a vi- 
able description for multiple state transitions of a black hole binary without invoking 
radial motion of the inner accretion disc. 

Key words: accretion, accretion discs - black hole physics - X-rays: binaries - X-rays: 
individual (GX 339-4) 
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1 INTRODUCTION 

The nature of the accretion flow for black hole binaries 
(BHBs) undergoing state transitions remains unclear and is 
an active topic of observational and theoretical pursuit. The 
phenomenology describing the predominant spectral states 
experienced by BHBs include the high/soft state (HSS), 
characterised by high luminosity and a standard thin accre- 
tion disc as the dominant spectral component; the low/hard 
state (LHS), characterised by low luminosity with a domi- 
nant hard power-law spectral component that is likely pro- 
duced by inverse Compton scattering of soft p hotons by a 
relativistic electron corona (|Shapiro et al.lfl976l ) : and the in- 
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termediate state (IS), which is itself further divided into the 
hard intermedi ate state (HI S) and soft intermediate state 
(SIS) (see, e.g.. lBelloni|[20T0l for an extensive description of 
these phenomenological states). 

Currently, the picture of a black hole accreting at a 
significant fraction of the Eddington mass accretion rate, 
M^dd, from a geometrically thin, optically thick disc ex- 
tending to the innermost stable circular orbit (ISCO) is the 
accepted description of the HSS. However, there is a lack 
of consensus regarding a universal geometry of the disc and 
innermost accretion flow in the LHS when the black hole 
feeds at a low fraction of A/Edd and/or the flow has a low 
ra diative effic i ency. According to the paradigm put forth 
bv lEsin et all j 19971 ), the accretion disc in the LHS is radi- 
ally recessed (i.e., truncated) with the interior region being 
composed of a radia tively inefficient accretion flow (RIAF) 
)|Naravan fc Yi| [l994). Alternatively, models where the inner 
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disc remains at the ISCO in the LHS have been proposed by 
invoking mag netically reconnect ing flares driving a coronal 
outflow (e.g.. lBeloborodovlll999h . Another model advocat- 
ing against disc truncation is a standard thin disc that gen- 
erates a magnetically dominated corona where the amount 
of accretion energy liberated in the corona increases wit h 
decreasin g mass acc r etion rate (jMerloni fc FabianI |2002h . 
Recently, Reis et al] |20l3) proposed that the intermedi- 
ate states may differ from one another only in the geomet- 
rical extent of the X-ray corona, with the SIS being the 
result of a collapse of the previously extended corona as 
the system moves into the disc-dominated HSS. The notion 
of an accretion disc, where the magnetorotational instabil- 
ity governs dissipation and angular momentum transport 
(Balb uTfe Hawlevl [l991). producing a magnetically domi- 
nated corona is supported by sophi sticated numerical simu - 
lations of local sheari ng boxes (e.g.jHirose et al.l l2006. 2009) 
and global discs fe.g.. lBeckwith et al.l2009l h By varying only 
the mass accretion rate, the spectral signatures of the HSS 
and LHS were reproduc ed in the global disc simulations of 
ISchnittman et all ([2012) , without the need for radial vari- 
ation of the inner accretion disc. The physical mechanism 
responsible for triggering a state transition is poorly under- 
stood, making it difficult to conceive what would drive a 
sudden disc recession or a drastic change in coronal activ- 
ity required to explain the observed spectral evolution; thus, 
lending credence to the need for a focused study of disc evo- 
lution during state transitions. 

Observational studies of BHB spectra have success- 
fully applied a multi - colour disc (MCD ) blac kbody model 
(|Mitsuda et al.l Il984l ; iMakishima et al.l Il986t ) to the ob- 
served quasi-thermal soft X-ray spectral comp onent in the 
HSS and even in the LHS in some cases (e.g., iMiller et al.l 
120061 : iReis et alJlioioh . Provided the observational param- 
eters are well-constrained, the radial location of the inner 
disc boundary can be measured from the best-fit MCD 
model parameters. However, this technique makes the crit- 
ical assumption that the deviations to the emergent disc 
spectrum caused by the disc vertical structure can be 
parametrized by a frequency-independent constant, known 
as the colour- correct ion, or spectral hardening, factor, / co i 
(|Shimura fc Takaharalll995h . 

Recently, investigations into the possibility of a vari- 
able /coi experienced a strong revival in the observational 
community. Studies of black hole X-ray binary accretion 
discs across various luminosity states using d ata from the 
Ross i X-ray Timing Explorer (RXTE) (e.g., Nowak et al 
20021; llbragimov et al.ll2005l ; IWilms et aLll2009 ; IDunn et al 



201 ll ) and the Swift observatory fRcvnold s~fc Miller! l201lT ) 



provide extensive systematic analyses of the properties of 
disc e voluti on in BHBs. The g l obal studies of IDunn et al.l 
(|201ll ) and IRevnolds fc Miller! (|201lT ) examined all of the 
available archived data for all of the well-known, confirmed 
BHBs. For nearly all BHBs where an outburst was observed, 
IDunn et alj (|201ll ) found that / co i remained relatively con- 
stant in the most disc-dominated states (i.e., the HSS), 
closely following the theoretical Ldi sc oc relation between 
the disc luminosity and disc effective temperature. For obser- 
vations with a relatively weaker disc component, character- 
istic of the IS, deviations from the luminosity-temperature 
relationship arose. Supposing that the inner disc radius does 
not change and allowing for a variable colour correction fac- 



tor, 1.6 < /coi < 2.6, where / co i increased as the disc frac- 
tion decreased, the Ldisc oc scaling was recovered. This 
finding sets the stage for a thorough consideration of an in- 
creasing /coi when a BHB transitions from the HSS to the 
LHS as a possible alternative to a truncat ed disc geometry in 
the L HS. Performing a similar analysis, IRevnolds fc Miller] 
(|201ll ) speculate that / co i may attain values as high as 5 as 
the BHB spectrum becomes less disc dominated. The main 
differences between our work and these recent observational 
studies are: (1) We focus exclusively on a single source, the 
well-known BHB transient, GX 339-4, over multiple state 
transitions. (2) We systematically explore different Comp- 
tonisation models (3) We introduce a novel approach for 
making relative measurements in order to reduce uncertain- 
ties. 

The methodology of this paper is to systematically in- 
vestigate to what extent variations in the disc structure 
and/or disc truncation can explain the observed spectral 
evolution for multiple state transitions of a BHB. Most im- 
portantly, we will demonstrate that a colour correction fac- 
tor that dynamically evolves within a physically reasonable 
range provides an adequate alternative to a truncated disc 
for explaining the changing disc spectrum. 



2 OBSERVATIONS AND DATA REDUCTION 

Motivated by the lack of consensus on accretion disc evolu- 
tion over the course of a BHB state transition, we aim to 
address the question: Can physically reasonable changes in 
the vertical disc structure alone explain the observed evo- 
lution of disc properties during a state transition without 
invoking a truncated disc geometry? In light of this ques- 
tion, we sought a source that met the following criteria: 
(1) The source was observed over multiple state transitions, 
permitting comparison of disc variations between different 
transitions. (2) Given that our methodology hinges on the 
presence of a disc, a disc component must have been ro- 
bustly detected in the LHS. (3) Conflicting claims for both 
a severely recessed disc and a disc extending to the ISCO 
in the LHS are published for the source. (4) Variable low- 
frequency quasi-periodic oscillations (LFQPOs), which may 
be linked to the disc evolution, have been observed coinci- 
dent with state transitions. (5) Reasonable constraints are 
known on mass, distance, and inclination; however, this is 
of ancillary importance. Based on these criteria, we selected 
GX 339-4 as the focus of our analysis. 



2.1 GX 339-4 



GX 339-4, discovered by iMarkert et al.l (|l973t ). is a low- 
mass X-ray binary and recurring transient hosting a dy 
namically confirme d black hole with a mass M > 5 .8Mq 
(|Hvnes et al. I l2003l : see also 
distance D = 8 ± 2 kpc 



Munoz-Darias et al.ll2008l ) and 
Zdziarski et al. 20041: see also 



iHvnes et al.l 12004^ . Supposing that the inner disc inclina- 
tion, i (0° for face on, 90° for edge on), can be assumed 



1 Reynolds & Miller (2011) also consider both phenomenological 
and physical Comptonisation models. 
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Trans. ID 


Start 
dd/mm/yy (MJD) 


End 

dd/mm/yy (MJD) 


Type 


JVobs 


^y-pow/ bknpow 
disc 


disc 


disc 


,,min 
QPO 

(Hz) 


, .max 
QPO 

CHzl 


"Kcp ' l Kep 


R02 


18/04/02 l'52SS2'i 


23/05/02 C524171 


Rise 


33 


32 


22 


31 


159fl1 a 


7 9(4 ) a 


14 


D03 


15/02/03 (52685) 


11/05/03 (52770) 


Decay 


53 


33 


30 


23 


0.29(l) a 


8.8(3) a 


9.7 


R04 


05/07/04 (53191) 


25/08/04 (53242) 


Rise 


30 


21 


21 


21 


0.306(7) a 


8.0(2) a 


8.8 


R07 


17/01/07 (54117) 


23/02/07 (54154) 


Rise 


34 


28 


17 


29 


0.142(3) a 


6.67(2) a 


13 


D07 


02/05/07 (54222) 


06/06/07 (54257) 


Decay 


46 


10 


39 


40 


0.070(3) a 


4.05(9) a 


15 


RIO 


26/03/10 (55281) 


20/05/10 (55336) 


Rise 


59 


11 


9 


10 


0.22(1) 6 


6.7(2) b 


9.8 


Dll 


07/01/11 (55568) 


15/02/11 (55607) 


Decay 


28 


12 


12 


12 


0.9(l) c 


4.5(3) b 


2.9 



Table 1. Table summarising the key features of each GX 339-4 state transition. From left to right, the columns are the state transition ID, 
start time of the transition, end time of the transition, type of transition, number of archived RXTE observations, number of observations 
requiring a disc component for Comptonisation models pow/bknpow, simpl, and comptt, minimum/maximum LFQPO frequency observed, 
and fractional change in inner disc radius implied if ^qp'q a nd Vq^q are associated with the Keplerian frequency and assumed to trace 
the inner disc location. ^Shaposhnikov fc Titarchukl j200Sh : iMotta et all 11201111 ; ^Stiele et al] d201lh 



to align with the binary inclina tion, then the lack of ob- 
served eclipses constrain i < 60° (ICowlev et alj|2002l). while 



radio observations imply i < 30° (|Gallo et al.l 2004). Re- 



sults from modelling the X-ray reflection features of GX 
339-4 in the S IS also suggest a l ower d isc inclination of 
i ~ 10° - 30° ([Miller et all 12004 l2009bl ) and in the LHS 
()Reis et al.ll200S ). Owing to the uncertainties on the intrin- 
sic parameters (namely, M, D, i) of GX 339-4 and the lim- 
itations of existing X-ray data, seeking absolute measure- 
ments that depend on these parameters presents a serious 
challenge. 

GX 339-4 experienced multiple outbursts since its dis- 
covery, making it an appealing source for temporal evolu- 
tion studies. The seven state transitions we select to study 
are listed in Table [T] The notation convention we adopt 
for referring to a transition is the letter R or D, signify- 
ing a rise or decay type, respectively, followed by the year 
associated with the transition. A rise (decay) type indi- 
cates a transition from the LHS (HSS) to the HSS (LHS). 
The intervals of each transition were selected based on pre- 
viously publis hed dates and the emergence/disappearance 
of LFQPOs JShaposhnikov fc Titarchukl 120091 ; IMotta et all 
120111 : IStiele et al.ll201lT )l 2 l We chose to extend the transition 
start/end times by ~ 5 days beyond the intervals quoted in 
the literature to ensure that we capture the full transition 
behavior in our subsequent analysis. 

Alarming discrepancies in the location of the inner ac- 
cretion disc edge, particularly in the LHS, appear in the 
literature based on spectral fits to X-ray observations of GX 
33 9-4. The strongest c ase for disc truncation is presented 
in iTomsick et al.l (|2009h , which used fluorescent iron emis- 
sion lines to find R[ n > 175 R g (90% confidence, assuming 
i = 30°) when GX 339-4 was in the LHS with luminosity 
Lx — 1.4 x 10 -3 Z/Edd (assuming D — 8 kpc). Claims of a 
truncated disc in the LHS also appear based on modelling of 
direc t disc emission accou nting for irradiati on of the inner 
disc (Cabanac et al.ll2009l ). Even in the SIS. lYamada et ail 
(|2009F suggest a recessed disc with Ri n — 5 — 32 R g (68% 
confidence, assuming i = 30°). These measurements sug- 
gest that the inner portion of the disc is absent in low 
luminosity states, perhaps replaced by a RIAF. In stark 



2 The end date listed for transitio n ID GX339-D03 in Table 1 of 
IShaposhnikov fc Titarc huk (200j|) is mistyped and should read 
06/05/03 (June 05 2003). 



contrast, there are numerous claims of the inner disc ra- 
dius in GX 339-4 remaining consistent with the ISCO in 
the LHS. Using physically motivated, sophisticated spec- 
tral fitting techniques, the disc inner radius wa s determined 
to lie within R in = 2.0 - 3.0 R g in the SIS (|Miller et al.l 
12004 ) and R in = 3.0 - 5.0 R g in the LHS w ith luminos- 
ity L x — 0.05 Z/Edd (assuming D = 8 kpc) (|Miller et al.l 
2006). Both of these results were confirmed by a systematic 
reanalysis that e mployed detailed disc reflecti on modelling 
(|Reis et al. 200 cj). Using the same data set as in lMiller et al.l 
(|2006T l, I Wilkinson fc Uttlevl (|2009l ) also measured a disc ex- 
tending down to 4 R g from their best-fit iron line model. 
ITomsick et"aH l|200Sf ) observed GX 339-4 in the LHS at lu- 
minosities Lx — 0.023 Z/Edd and Lx — 8 x 10 -3 Lsdd, find- 
ing the disc to extend within 10 R g an d being consistent 
with the disc remaining fixed at ~ 4 R s . iReis et al.l (|2010h 
performed a systematic study of eight BHBs (including GX 
339-4) in the LHS, robustly detecting a thermal component 
in all sources and placing a stringent upper limit of 10 R g 
on the truncation radius for six sources confirmed by both 
broad iron line and thermal disc modelling independently. 

In order to measure black hole spin, one must associate 
Rin with the ISCO, located at Ziisco- A measurement of 
the inner disc radius is the crucial diagnostic for measuring 
black hole spin, as Ziisco is a monotonically decreasing func- 
tion of the dimensionless spin parameter, a * = Jc/GM , 
whe re J is the black hole angular momentum (|Bardeen et al.l 
1972). Perhaps unsurprisingly, conflicting claims regarding 
the behaviour of Ri n have led to discrepancies in spin mea- 
surements of GX 339-4. Spin measurements e mploying the 
relativistically b roadened iron line method (|Fabian et al.l 
119891 : lLaorl Il99ll ) suggest that GX339-4 harbours a near- 
maximally spinning black hole with a, = 0.94 ± 0.02 
(|Miller et al.l l2009bl ; IReis et all 120081), while the th ermal 
disc continuum fitting technique (Z hang et al.l [l997) sup- 
orts a low-to-moder ately spin ning black hole with a, < 0.4 
Yamada et all 120091 ; see also iKolehmainen k, Done! |2010| ; 
iKolehmainen et"aL 201ll ). 

LFQPOs exhibit variability si multaneous with trans i- 
tions between the LHS and HSS |Rutledge et all Il999bl ). 
Variable LFQPOs were observed in each GX 339-4 transi- 
tion studied here, with the measured frequenc y ranges listed 
in Table [J and typically spanning 0.2 — 8 Hz jBelloni et al.l 



20051; IShaposhnikov fc Titarchukl |2009|; IMotta et al.l l201ll : 



Stiele et al.ll201ll ; lNandi et al.ll2012t ). The idea that LFQPOs 



may originate from oscillating modes of a quasi-spherical 
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RIAF that is interior to, or partially overlapping with, 
a trun cated thin disc was proposed by iGiannios fc Spruit) 
(2004). Long time-scale variability of the disc blackbody 
component of GX 339-4 was associated with disc instabili- 
ties and interpreted as arising f rom intrinsic variability of the 
disc ijWilkinson fc Uttlevi 2009). This suggests that the vari- 
able LFQPOs may be intricately linked to disc truncation; 
however, the truncation radius cannot exceed ~ 20 R g if the 
observed variations are limited by the viscous timescale. 

2.2 Rossi X-ray Timing Explorer 

The monitoring campaign of the Rossi X-ray T iming Ex- 
plorer (RXTE) l|Bradt et al.l Il993l ; ISwankl Il999l ^ observed 
GX 339-4 transitioning between states in 2002, 2003, 2004, 
2007, 2010, and 2011, obtaining well-sampled (in time) ob- 
servations spanning each transition before being decommis- 
sioned in early 2012. We examined 283 RXTE observations 
of GX 339-4 covering seven time intervals, during which 
the source experienced state transitions (see Table [TJ. The 
RXTE standard products, comprised of the source and back- 
ground spectra along with instrument response matrices, 
were obtained from the High Energy Astrophysics Science 
Archive Research Cent er data archive__ T he Proportional 
Counter Array (PCA) dJahoda et al.lll996l) and High En- 
ergy X-ray Timing Experiment CHEXTE) |Rothschild et al.l 
1998) standard source spectra were analysed jointly. The 
PCA energy spectrum (Standard 2 mode) is built from sum- 
ming all layers from various Proportional Counter Units 
(PCUs). The PCU-1, PCU-2, and PCU-3 are used in all ob- 
servations prior to 25 December 2006 (i.e., R02, D03, R04), 
while only the PCU-2 and PCU-3 are used in observations 
thereafter (i.e., R07, D07, R10, Dll) due to a propane loss 
in the top layer of PCU-lQ For observations corresponding 
to transitions R02, D03, R04, R07, and D07, we used the 
HEXTE Cluster B spectrum, but could not use these data 
for R10 and Dll because the HEXTE Cluster B was perma- 
nently left in an off-source pointing position from 29 March 
2010 onwardQ The HEXTE Cluster A data for transitions 
R10 and Dll are affected by strong residuals, resulting from 
difficultly in determining the appropriate contribution of the 
background spectra. Consequently, we excluded all HEXTE 
data while analysing transitions R10 and Dll, using PCA 
data only for these observations. This resulted in many ob- 
servations with a strong soft component, which is indicative 
of a disc, being excised from our sample due to poor spectral 
modelling of the high-energy continuum (see t]4.1[) . 

While the time sampling capability of the RXTE satel- 
lite makes it well-suited for a study of BHB state transitions, 
which are a highly time-dependent phenomenon, RXTE has 
important limitations. Determining the spectral properties 
of the disc component for a given observation is made dif- 
ficult by the limited low-energy coverage of RXTE, which 
only extends down to ~ 3 keV, while the disc spectrum 
peaks at ~ 1 keV. Especially in the HIS and LHS, the dimin- 
ished disc component, if present, is difficult to characterise 
due to a combination of its weakness and the dominance of 

3 http:/ /heasarc. gsfc.nasa.gov/cgi-bin/W3Browse/w3browse. pi 

4 http:/ /heasarc. nasa.gov/docs/xte/recipes/stdprod_guidc. html 

5 http:/ /heasarc. nasa.gov/docs/xte/xhp_new. html 
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Figure 1. Hardness-intensity diagram for all observations of GX 
339—4 considered in this work generated from the results of the 
pow/bkn model (see i|4.2H . The disc fraction for each observation is 
given by the glyph colour and is an indicator of the spectral state. 
Data points marked by a grey cross correspond to observations 
that were not deemed to statistically require a disc component. 
The key in the bottom left identifies the glyph associated with 
each transition ID . See i|2.2l for a description of how the HID was 
generated. A wide range of hardness is considered in this work. 



a Comptonised component. Figure [T] presents a hardness- 
intensity diagram (HID) of all 283 RXTE observations of 
GX 339-4 considered in this work. Observations that do not 
statistically require a disc component (see > )4.1I) are indi- 
cated by grey crosses, while the various glyphs indicate the 
transition associated with the observation. The plethora of 
IS and HSS observations (i.e., low-to-moderate hardnesses) 
deemed to not require a disc comes from the RIO and Dll 
data sets, where the absence of HEXTE data prevented ac- 
ceptable fits to disc-dominated spectra. Notably, Figure [T] 
shows that the set of observations considered in the subse- 
quent analysis span all spectral states, which indicates that 
we are studying the complete evolution of state transitions 
in GX 339-4. 

For a given observation, hardness and intensity are cal- 
culated from the PCA data as follows. The distribution of 
counts across all channels is obtained from both the stan- 
dard PCA source spectral file and background spectral file. 
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The background counts are then subtracted from the source 
counts channel-by-channel, yielding the counts in each chan- 
nel associated with GX 339-4. The hardness is then defined 
as the ratio of GX 339-4 counts in channels 15 - 24 (6.3 
- 10.5 keV) to tho se in channels 8 - 14 (3.8 - 6.3 keV) 
(|Belloni et al.l 120051 ). A count rate is obtained by summing 
the counts across all channels and dividing by the exposure 
time. Good time intervals (GTls) are identified from the GTI 
file included in the standard products. The average number 
of PCUs that are turned on during the GTIs is determined 
from the filter file included in the standard products. The 
intensity is then defined as the count rate divided by the av- 
erage number of active PCUs during the observation GTls. 

To justify using the \ 2 fitting statistic for deriving pa- 
rameter constraints, both the PCA and HEXTE spectra 
were binned to require a minimum of 20 counts per en- 
ergy bin using the FT00L grppha. A systematic error of 0.6% 
was added to the PCA spectrum using grppha; however, the 
PCA spectra without systematic errors added are used when 
deriving relative measurements (see £|4.3.1|). All spectr al fits 
are performed with XSPEC version 12.7.1 <|Arnaudll996h . The 
PCA and HEXTE spectra are fitted over the 2.8 - 25 keV 
and 20 — 100 keV energy ranges, respectively. The response 
energy range over which spectral models are calculated is ex- 
tended to 0.1 keV and 100 keV using the energies extend 
command in XSPEC with 50 logarithmic steps in each direc- 
tion. All observation times correspond to the midpoint of 
the PCA observation. Quoted errors from spectral fits are 
la uncertainties, unless specified otherwise. 



3 MULTI-COLOUR DISC BLACKBODY 
MODELS 

The accretion disc in BHBs manifests itself as a soft, quasi- 
blackbody component w ith a typical peak inner disc temper- 
ature of ~ 1 keV (e.g., iFender fc Bellonil 120041 ). While the 
accepted description for the origin of the soft thermal com- 
ponent in the HSS is a geometrically thi n, optically thick 
accretion disc (|Shakura fc Sunvaevlll973l) extending down 
to the ISCO, no such consensus has been reached for the 
nature of the innermost flow in the HIS or LHS. Appealing 
to spectral fits of the thermal continuum with multi-colour 
disc (MCD) blackbody models to isolate the disc component, 
one can measure the temporal evolution of the accretion disc 
over the course of a state transition. 



3.1 The Standard MCD Model 

The most commonly adopted disc model is the MCD 
blackbody model ( dis kbb in XSPEC) (|Mitsuda et al.l Il984l ; 
iMakishima et al.l Il986h . The attractiveness of the MCD 
model is owed to its legacy, applicability, and simplicity. 
However, the familiarity and broad acceptance of diskbb 
often mask subtle features of the model, possibly leading to 
unintended misuse. Here, we dissect the diskbb model in 
order to elucidate its underlying assumptions and expose its 
limitations. Assuming a geometrically thin, optically thick 
accretion disc that locally emits a blackbody spectrum in 
concentric annuli, the colour temperature for a MCD model 



at a given radius is described as, 

-3/4 



T col (r) = T, 



1-/3 



Rin 



-1/2 



1/4 



(1) 



where R[ n is the radius of the inner disc edge and /3 describes 
the boundary condition at this location. For discs in the 
HSS, -Ri n is commonly associated with Risco- A zero-torque 
boundary condition at R[ n corresponds to /3 = 1, while /3 = 
allows for non-zero stresses at the inner disc and is the choice 
adopted in diskbb (|Zhang et al II 19971 ; iGierliriskT et al II 19991 ; 
IZimmerman et al.ll2005l ). The leading factor, which roughly 
describes the colour temperature at the inner disc, is given 
by, 



T* = U 



f 3GMM \ 

\~^rl) 



1/4 



(2) 



where G is the gravitational constant, a is the Stefan- 
Boltzmann constant, M is the compact object mass, M is 
the m ass accretion rate, and ,f co i is the colour correction 
factor ( Sh imura fc Takaharai ri995). also referred to as the 
spectral hardening factor. For a zero-torque condition at Ri n 
(i.e., j3 = 1 in Equation [1]), T» is related to the maximum 
temperature attained in the disk by Ted, max = 0.488T,, 
which occurs at r = (49/36) i? m . Since the disc temperature 
does not reach T, at any radius, Ted, max is a more physi- 
cally meaningful characterisation of the blackbody temper- 
ature associated with the spectrum. For a thin disc that is 
torqued at the inner boundary, as is the case for diskbb, 
Tco^max coincides with T„ because the temperature pro- 
fil e grows asymptotically with decreasing r (see Figure 1 
of lZimmerman et al.ll2005h . One of the two free parameters 
of the diskbb model is T, , implicitly assuming / co i = 1. The 
'normalisation' is the other free parameter of diskbb and is 
given by, 



100 



/col 



Rill 



D 



kpc 



(3) 



where, again, / co i = 1, Dk P c is the distance to the source in 
kiloparsecs, -Rm.km is the inner disc radius in kilometers, and 
i is the inner disc inclination. Inner disc radii are commonly 
derived from Khb by measuring intrinsic system properties 
and imposing a constant / co i (not necessarily unity as as- 
sume d by diskbb), regardle ss of the source luminosity state 
(e.g., IZdziarski et al.l 120041 ). Remarkably, the simple two- 
parameter diskbb model fully describes the accretion disc 
spectrum under the appropriate assumptions of a non-zero 
torque inner disc boundary condition and a fixed colour cor- 
rection factor for characterising the overall deviation from a 
blackbody spectrum. 



3.2 Zero- Torque Boundary Condition and Disc 
Vertical Structure 

If the inner disc radius can be identified with the ISCO, 
a zero-torque boundary condition becomes attractive from 
the physical consideration that the gas viscous inspi- 
ral timescale is long relative to the free-fall timescale 
for gas at the ISCO; therefore, the gas interior to R[ n 
does not ha ve time to radiate appreciably before be- 
ing a c creted |Shakura fc Sunvaevlll973l; iNovikov fc Thornel 
Il973l ; IPage fc Thornel 1 19741 ; lAbramowicz fc Katol 1 19891 ; 
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iFrank et al.l Il992l ; lAfshordi fc Paczvnskil 120031 ). In princi- 
ple, for a geometrically thin disc in the presence of a 
magnetic field, non-negligible torque resultin g from mag- 
netic stresses may be present at the ISCO (Krolik 1999; 
Gammiel 1 19991 : lAgol fc Kroiikl |200d : lAfshordi fc Paczvnskil 



20031 ). However, magnetohydrodynamic three-dimensional 



simulations of global thin accretion discs around black holes 
demonstrate that a weak magnetic field is unable to couple 
the gas within the plunging region to the disc, despite large 
fluctuating magnetic stresses inside the ISC O, in accordance 
with the zero-torque boundary c ondition l|Armitage et al.l 
l200ll : iRevnolds fc Armitagel l200ll ). Similar simulations ex- 
ploring initially purely poloidal and toroidal field topologies 
instead observe significant stress at the ISCO and continu- 
ing deep wit hin the plunging region, refuti ng the zero-torque 
assumption jHawlev fc Krolikll200ll . l2002l ). A difficulty with 
allowing for a non-zero torque in disc models is that there 
is no obvious choice for the magnitude of the torque, ne- 
cessitating the introduction of a free parameter. The zero- 
torque choice is both physically motivated and is not ar- 
bitrary; however, the presence of significant torque on the 
inner disc may be important in real systems depending on 
highly unknown factors such as the disc thickness, dissipa- 
tion properties, and magnetisation. 

The underlying physics responsible for altering the lo- 
cally emitted spectrum of an optically thick disc from a pure 
blackbody are parameterised into the colour correction fac- 
tor. The result of this approximation is the colour-corrected, 
or 'diluted', blackbody, 



/col 



Bv (/colTeff) 



(4) 



where I v is the specific intensity, B„ is the Planck func- 
tion, T c ff = T co i//coi is the effective temperature, and 
/coi = 1 refers to a pure blackbody spectrum. In re- 
ality, deviations in the emergent disc spectrum from a 
canonical blackbody arise from frequency-dependent opac- 
ities determined from the disc vertical structure, which 
is neglected in the standard thin disc treatment. The 
degree of spectral modification will be governed by the 
combined effects of electron scattering (i.e., Compton- 
isation in the disc), free- free emission/absorption, and 
bound-free absorption, resulting in a departure from lo- 
cal therm odynamic equilibrium and a hardening of the 
spectrum dFelten fc Reeslll972l; IShimura fc Takaharalll995l: 



Zavlin et al.lll996l : iRutledge et al.lll999al : iMcClintock et al.l 
20041 ). A surrounding corona acting as a depository for 



some fraction of the disc accretion power provides yet 
another means for altering the emergent disc spectrum 
(|Svensson fc Zdziarskil [l994l : iMerloni et al.ll2000l ). Further- 
more, magnetic pressure support likely contributes signifi- 
cantly to hydrostatic balance in the disc photosphere and 
may act to vertically extend the disc a t mosphere, producing 
a har der spectrum ( Bla es et al.l 120061 : iBegelman fc Pringld 
2007). The phenomenological colour-correction prescription 
avoids all of these complications by neglecting any frequency 
dependence, only permitting translational (i.e., frequency- 
independent) hardening and attenuation of the Planck func- 
tion. Provided that this model yields a satisfactory fit to the 
observed disc spectrum, one presumes that / co i , which is as- 
sumed to be constant for all radial locations in the disc, 



provides an adequate description of the disc vertical struc- 
ture in an average sense. 

The two disc models we appeal to in this work dif- 
fer from diskbb in that they both adopt a zero-torque 
inner boundary condition and allow for the possibility of 
a variable colour correc tion factor. The ezdiskbb model 
(|Zimmerman et al"1 [2005) is the analog of diskbb, but with 
the zero-torque boundary condition enforced. The free pa- 
rameters are the maximum disc colour temperature, T co i imax 
defined above, and the normalisation, given by 



100^ 
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■Riii 



/col V ^kpc 



(5) 



where / co i is not specified. We emphasize that changes in 
.Kez do not necessarily imply changes in the disc inner ra- 
dius, due to the param eter degeneracy betw een Ri n and 
/ co i. The diskpn model (Gicrlinsk i et al.lll999l ). along with 
the zero-torque boun dary condition, incorporat es a pseudo- 
Newtonian potential (|Paczvriskv fc Wiitalll980l ) to more ac- 
curately calculate the radial temperature distribution in the 
disc at the expense of introducing the inner disc radius as 
a free parameter. The two remaining free parameters are 
Tcoi.max and the normalisation, 
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/col 
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where Mm is the compact object mass in solar mass units. 
In all of the MCD model normalisations, / co i enters to the 
fourth power. Therefore, asking how much f co \ could con- 
ceivably change during a state transition is not unreason- 
able, as only small changes in / CD i are required to explain 
the observed spectral evolution, which may not necessitate 
disc truncation. 

The disc models described up to this point neglect 
many important physical considerations. We wish to em- 
phasize that sophisticated, physically motivated models are 
available in XSPEC that accoun t for irradiation of the disc 
(diskir; iGierlihski et al.ll200Sl ). general relativistic e ffects 
around a Schwarzschild (grad; 



Kerr black hole (kerrbb; iLi et al 



Ebisawa et al.1 fl991) and 
20051) . and detailed mod- 



elling of t he disc vertical st ructure around a Kerr black hole 
(bhspec; iDavis et al.ll2005l ). Realistic disc models come with 
the expense of introducing additional parameters. The lim- 
ited quality of X-ray observations and the desire to limit the 
number of model parameters forces the use of over-simplified 
disc models, with the hope that they can adequately cap- 
ture the most important aspects of the accretion physics. 
The purpose of the subsequent analysis is to demonstrate 
that simplistic disc models produce good spectral fits but 
they cannot be used to make reliable claims regarding the 
behaviour of the inner disc radius. 



4 ANALYSIS AND RESULTS 

In what follows, we exploit the timing and monitoring capa- 
bilities of RXTE, which provide multiple sequential spectra 
of moderate resolution tracking the transition between spec- 
tral states, to study the evolution of the accretion disc in GX 
339-4. The ultimate goal of our spectral analysis is to iso- 
late the disc component in order to measure relative changes 
in disc properties during transitions over the rise (LHS — >• 
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HSS) and decay (HSS — ► LHS) stages of an outburst. We 
show that allowing for physically reasonable changes in f co \ 
provides an alternative explanation to disc truncation for 
the observed disc spectral evolution during state transitions. 
Table [3] provides a sample of best-fit parameters for the fits 
used in this work. Table in its entirety is available as an 
electronic supplement. 

4.1 Requirement of a Disc Component 

We wish to verify the presence of a thermal disc compo- 
nent with adequate confidence using phenomenological mod- 
els that have proven successful in characterising the gen- 
eral features of BHB X-ray spectra. To determine whether 
or not a disc component is required in a given RXTE 
spectrum, we adopt a methodology based on x 2 statistics 
that uses phenomenological models to fit numerous obser- 
vations in various spectral states , as motivated by previ- 
ous w orks |Dunn et al.ll2008l . l2010l . l201ll ; lRevnolds fe Milled 
120111 1. The PC A and HEXTE data are fit jointly, allow- 
ing an energy-independent constant factor to vary for the 
HEXTE data to account for any normalisation offset rel- 
ative to the PCA and linking all parameters between the 
PCA and HEXTE spectra. The spectral models are mod- 
ified by photoelectric absorption using the phabs model 
with a fixed equivalent neutral hydrogen co lumn density of 
N H = 5.7 x 10 21 cm -2 jMiller et al.ll2009tj ). Both Galactic 
and intrinsic absorption are folded into the phabs model, 
leading to a potential ambiguity as to whether low-energy 
spectral evolution is due to chan ges in the intrinsic a bsorp- 
tion or to evolution of the source. iMiller et al.l (|2009al ) fitted 
high-resolution BHB spectra (including GX 339-4) over a 
wide range in luminosity and spectral states, demonstrat- 
ing that the low-energy spectral evolution is a consequence 
of the evolution of the source itself, rather than a variable 
neutra l photoelectric absorption. Notably, ICabanac et al.l 
(2009) present the case for possible variations in Nu during 
outbursts of GX 339-4; however, the restricted low-energy 
coverage of RXTE does not permit Nh to be a free param- 
eter in our spectral fits. 

As a first attempt, the spectra are fitted with a simple 
model consisting of an absorbed power-law and MCD black- 
body (i.e., phabs x (powerlaw+diskbb) in XSPEC) to model 
the continuum and accretion disc, respectively. In many ob- 
servations, this simple prescription does not yield an ade- 
quate fit due to the presence of a broad Fe Ka emission 
line. To check for the requirement of an iron line com- 
ponent, a Gaussian line profile (gauss in XSPEC) is added 
to the model at fixed centroid energy of 6.4 keV with the 
width, cr, and normalisation, Ku nc , as free parameters. An 
F-test is inappropriate for determining whether the inclu- 
sion of a spectral line is a st atistically significant improve- 
ment JProtassov et al.| 20021). Th e alternative diagnostics 
outlined bv lProtassov et al.l (|2002T l for testing for the pres- 
ence of a line are not easily incorporated into the spectral 
fitting framework used here. We are concerned with obtain- 
ing an acceptable fit in order to determine the presence of a 
disc component and do not wish to make any strong state- 
ments regarding Fe Ka e mission; ther e fore, we adopt the 
line detection procedure of iDunn et all l|200Sf ). In order for 
a line to be considered significant, two criteria must be met: 
(1) An F-test must determine that the line is significant at 



> 3<r (F-statistic probability V < 0.0027). (2) Given the la 
uncertainty in the line normalisation, <7K linc , obtained with 
the XSPEC error command, we require Kn ne > 3<JK line ■ This 
procedure for line determination is more conservative than 
appealing to an F-test alone, though we note that some de- 
tected line features may be erroneous. 

We next test for the presence of a high energy break in 
the spectrum by replacing the power-law component with 
the broken power-law model, bknpower, and repeating the 
aforementioned fitting procedure. If an F-test determines 
that the broken power-law model is a significant improve- 
ment over the power-law model at > 3<r confidence, then we 
proceed with the broken power-law model imposing the pa- 
rameter ranges for the spectral indices and break energy of 
1 < r < 4 and E^ > 10 keV, respectively. In instances where 
a line component was deemed necessary in the power-law 
model but not the broken power-law model, or vice versa, 
we select the best fitting model determined by the lowest 
reduced x 2 -statistic, xl = X 2 ' l v i where v is the number of 
degrees of freedom in the model. For cases when the best-fit 
converges to the parameter limits of the bknpower model, the 
pipeline reverts to the powerlaw model for that observation. 
If the best-fit powerlaw model also runs into a parameter 
bound, the observation is discarded. 

At this point in the phenomenological fitting pipeline, 
for a fit to be considered acceptable for use in subsequent 
analysis, we require xi < 2. Often, a poor fit to the high- 
energy spectrum is the culprit for failure to meet this crite- 
rion. For instance, a poor fit may result from the inability 
of simplistic, unphysical models to properly fit a reflection 
component or a complicated spectrum arising from a dy- 
namic corona. Since we are using phenomenological models 
to characterize the soft disc component and are not inter- 
ested in a detailed characterization of the high-energy spec- 
trum, if the best-fit model to the combined PCA/HEXTE 
data gives Xm > 2, then we repeat the above fitting proce- 
dure on the PCA data alone in search of an acceptable fit. In 
these cases where only the PCA data are being fitted, we use 
the powerlaw model only because any spectral energy break 
is expected to lie beyond the high-energy cutoff of the PCA, 
rendering the use of the bknpower model inappropriate. 

The disc component (diskbb) is then removed from the 
best-fit model and the spectrum is re-fitted. An F-test is 
performed to determine whether or not the inclusion of the 
disc improved the fit by > 3<r. Observations that meet these 
criteria are kept for further analysis, while those that do 
not statistically require a disc are discarded. Running this 
model fitting procedure on 283 total observations of GX 339- 
4 in different stages of transition between the low/hard and 
high/soft states resulted in 177 observations where a disc 
component was deemed statistically required from spectral 
fitting techniques (see Table [1]). 

We caution that the absence of a quasi-blackbody com- 
ponent in a given RXTE spectrum does not imply the ab- 
sence or truncation of the physical accretion disc. The low 
energy cutoff of the RXTE PCA is restricted to ~ 3 keV; 
thus, only the Wien tail of the MCD blackbody spectrum 
will be in the observable band. Although we may not detect 
a disc, especially in the HIS and LHS, this is not enough 
to claim the absence of a disc component, which may be 
dominant at ene rgies below the sensitivity of RXTE (e.g., 
IDunn et al . 2008) and has been observed with XMM-Newton 
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in th e LHS of GX 339-4 (e.g.. iMiller et al.ll2006l ; iReis et al.l 
l2008h . 

Figure [2] shows the data-to-model ratios with the disc 
component removed from the best-fit model for representa- 
tive observations of each state transition. Disk fractions are 
defined by, 



^disc I ^total 



DF 



-.0.1-10 kcV 



7.0.1-100 kcV ' 



(7) 



where Fkisc and Ft ota i are the unabsorbed fluxes of the disc 
component and overall model, respectively, taken over the 
energy ranges indicated by the superscripts. Fluxes were ob- 
tained with the cf lux command in XSPEC. Disc components 
are robustly detected by the pipeline described above, even 
for moderately low disc fractions, which is essential for track- 
ing the evolution of the disc throughout each transition in 
the subsequent analysis. 



4.2 Sensitivity of Spectral Fits to the Continuum 
Model 

The motivation for employing a simple power-law model 
with limited free parameters is to fit the portion of the 
spectrum that is thought to arise from Comptonisation of 
thermal disc photons by a low density, energetic electron 
corona of unknown geometry, which produces a power-law 
spectrum to first order. The phenomenological power-law 
model extends to arbitrarily low energies, which does not 
mimic a physical Comptonisation process where a sharp cut- 
off is expected at an energy comparable to the characteristic 
energy of the input seed photon distribution, namely, the 
quasi-blackbody disc spectrum. Therefore, the extension of 
the power-law component toward arbitrarily low energies re- 
sults in 'flux stealing' from the disc component, making a 
disc more difficult to detect than had we adopted a more 
physical Comptonisation model. 

To investigate the sensitivity of fitting the disc compo- 
nent on the choice of Comptonisation model, we explore two 
different models to chara cterise the continuum. The Comp- 
tonis ation model comptt (|Titarchuklll994l ; lHua &: Titarchukl 
Il995t ) provides an analytic description of the spectrum pro- 
duced from a soft photon distribution propagating through a 
hot plasma cloud. We re-ran the disc detection pipeline with 
comptt in place of the power-law component, allowing only 
the plasma temperature, T p , plasma optical depth, r p , and 
normalisation, -Kcomptt, to be free parameters. The remain- 
ing comptt parameters that we chose were a redshift z = 0, 
a seed photon temperature, To = T co i jmax (i.e., 0.488 times 
the diskbb temperature), and a spherical plasma cloud ge- 
ometry. The additional free parameter introduced by comptt 
is justified by the physical, as opposed to phenomenological, 
nature of the model. For cases where the best-fit converges 
to the enforced plasma temperature limit, T p = 100 keV, 
the data are re-fitted with T p fixed at 100 keV. When the 
disc component is removed to check for the requirement of 
a disc, To becomes a free parameter, as there is no longer 
an obvious, physically motivated choice for the seed photon 
distribution injected into the corona in the absence of a disc. 
A total of 166/283 observations required a disc component 
when comptt was used, compared to the 177/283 using the 
standard phenomenological powerlaw and bknpower models. 
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Figure 2. Ratios of the RXTE data to the best-fit models with 
the diskbb component removed to highlight the need for a quasi- 
blackbody feature to model the soft excess. For each transition, 
three data/model ratios are shown that representatively span the 
range in measured disk fractions. Glyph colours denote the disc 
fraction for each fit. For comparison, the best-fit models where 
Comptonisation is modelled by the powerlaw or bknpower (left 
column) and comptt (right column) XSPEC models are shown. For 
each row, the observation IDs are given in the upper right corner 
of the plots and the transition ID on the rightmost ordinate. 
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Figure 3. Dependence of the strength and spectral shape of the 
disc component on the choice of Comptonisation model for tran- 
sition R07. The residuals plotted are computed by subtracting 
the data-to-model ratio for the pow/bknpow best-fit model from 
the best-fit comptt data-to-modcl ratio. The data-to-model ratios 
used are those where the disc component was removed, as done 
in Figure[2] Residuals falling on the dashed line indicate no differ- 
ence in the measured disk component between the two Compton- 
isation models, while a positive excess implies flux stealing from 
the disc component by the pow/bkn model relative to the comptt 
model. From top to bottom, the residuals were computed from 
observation IDs 92428-01-03-00 (red glyphs), 92428-01-04-02(Mue 
glyphs), and 92085-01-01-00 (green glyphs), which represent the 
HSS, IS, and LHS, respectively. The flux stealing effect is most 
pronounced in the HSS and IS. 

An alternati ve empirical subst itute for a power-law 
model is simpl l)Steiner et al.l [20091 ) . which truncates the 
power-law spectrum at an energy comparable to the char- 
acteristic input seed photon energy; thus, avoiding the un- 
physical low-energy extension of the powerlaw and bknpower 
models. Since simpl is a convolution model, we cannot ap- 
ply the disc detection pipeline of i|4.1l to determine whether 
or not a disc component is required. Instead, the previ- 
ously determined 177 best-fit models are re-fitted with simpl 
as the component responsible for modelling Comptonisa- 
tion, allowing for both up- and down-scattering of the input 
seed spectrum. Again, as was done for the disc detection 
pipelines using powerlaw/bknpower and comptt as the con- 
tinuum models, the observation is discarded if the best-fit 
converges on a parameter limit. This returned 150/177 ob- 
servations where the best-fit does not encounter a parame- 
ter limit. While simpl cannot capture the high energy break 
commonly observed in BHB spectra, often resulting in the 
dismissal of the HEXTE data in order to obtain acceptable 
fits, this excision of data is justified because a detailed model 
of the high-energy spectrum is ancillary to a reliable model 
of the disc spectrum. 

The choice of Comptonisation model influences the 
measured disc properties. As discussed above, the three 
Comptonisation models we explore in this work are: 
(1) standard power-law or broken power-law (hereafter, 
pow/bknpow), (2) power-law with a low-energy cut-off (here- 
after, simpl), and (3) physical Comptonisation model (here- 
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Disc Fraction (F disc / F total ) 

Figure 4. Histogram of disc fractions for fits performed with 
Comptonisation models pow/bknpow (grey filled), simpl (red for- 
ward slashes), and comptt (blue backward slashes). The excess 
of low disc fractions for the pow/bknpow model is a consequence 
of the power-law extending to arbitrarily low energies, effectively 
stealing flux from the disc component. 



after, comptt). For representative observations over all tran- 
sitions fit with the pow/bknpow and comptt models, Figure 
[2] shows data-to-model ratios with the disc component re- 
moved from the best-fit model, which is a probe of the spec- 
tral shape and strength of the disc component. The disc is 
strong in the HSS and gradually weakens, but remains de- 
tectable, in the LHS down to DF ~ 0.05. 

Figure [3] shows the residuals for the data-to-model ra- 
tios of the comptt fits compared to those for the best-fit 
pow/bknpow models for a representative transition. The flux 
stealing phenomenon inherent to the low-energy extension 
of the pow/bknpow model is apparent from the positive ex- 
cess in the residuals and affects the measured disc spectrum, 
particularly for high disc fractions that are characteristic of 
the HSS. A histogram of measured disk fractions for all ob- 
servations deemed to require a disk and for each choice of 
Comptonisation model is shown in Figure [4] Generally, re- 
placing the power-law model with more physically motivated 
Comptonisation models produces higher disc fractions and 
yields disc model spectra that are unplagued by flux stealing. 
The measured disk fraction for a given observation generally 
follows DF simp i 

Given this study on the sensitivity of the spectral fits to 
the choice of continuum model, we must decide how to pro- 
ceed with the subsequent analysis in order to best represent 
the true character of the disc evolution during a state transi- 
tion. In the interest of comparison and because discrepancies 
in the measured disc properties arise between the three dif- 
ferent, but seemingly adequate, Comptonisation models, we 
do not favour one model over another. Instead, we elect to 
perform the subsequent analysis in three branches, identical 
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Figure 5. Histogram of relative errors on the best-fit diskbb 
normalisation, o"A' bb /^"bbi f° r an RXTE observations requiring a 
disc component. Models fit with a systematic error of 0.6% added 
to all PCA energy channels and no systematic error added to the 
PCA data are shown as red forward slashes and grey bars, re- 
spectively. The relative error on the disc normalisation is reduced 
when systematic error is not added to the data. 



in all respects except for the Comptonisation model (i.e., 
pow/bknpow, simpl, comptt) being used in the spectral fits. 



4.3 Disc Evolution Results from Spectral Fitting 

A parameter degeneracy exists between the inner disc radius 
and the colour correction factor within the normalisation 
of the standard diskbb model. Therefore, decomposing the 
contributions to the disc evolution due to changes in the ver- 
tical disc structure, parametrized into / co i, and a migrating 
inner disc region is nontrivial. Combining the models diskpn 
and ezdiskbb, while appealing to relative measurements of 
the disc model normalisation evolution, we examine the de- 
gree of change in / co i and Ri u required to explain the disc 
spectral evolution in state transitions of GX 339-4. 



4-3.1 Relative Measurements 

Obtaining reliable absolute measurements of quantities 
locked up in the disc normalisation is limited by constraints 
on the observables, namely, the mass, distance, and disc in- 
clination, which all appear in the normalisation. However, 
relative changes in / co ] and/or 7?i n are accessible from the 
evolution of the disc normalisation because the intrinsic pa- 
rameters do not change between successive observations. 

Imposing a small systematic error on RXTE PCA data 
is common practice, resulting in a substantial improvement 
in the quality of spectral fits. Applying model fits to data 
where no systematic error is added will result in poor fits but 
with best-fit parameter errors dramatically reduced. While 



the best-fit parameter values may be unreliable, they should 
be affected in a systematic sense and systematic errors as- 
sociated with instrument calibration will be removed from 
the parameter errors. Provided that relative, rather than 
absolute, measurements are of interest, the evolution of a 
parameter may be tracked from fits to spectra without sys- 
tematic error incorporated into the PCA data. 

To test the validity of this proposition, we re-fitted 
all spectra deemed to require a disc component with the 
best-fit models determined from the aforementioned fitting 
pipeline, but with no systematic error added to the PCA 
data. Figure [S] compares the distribution of the diskbb nor- 
malisation relative error, defined as o-pc bb / Kbb, for fits to all 
GX 339-4 spectra when no systematic error and the stan- 
dard 0.6% systema tic error were added to the PCA data 
jjahoda et al.ll2006r i. As expected, not including systematic 
errors results in improved relative errors on the diskbb nor- 
malisation. Adopting an analysis technique based on omit- 
ting systematic errors from the PCA data requires that 
the disc normalisation values are only affected in a system- 
atic way as a result. For each transition, Figure [5] shows 
Kbb,rci = ifbb.nosys/Jfbb, which is defined as the best-fit 
diskbb normalisation for fits with no systematic error added, 
-Kbb.nosys, relative to the normalisation for fits with 0.6% 
systematic error added, Kbb- The constancy of this ratio 
throughout the spectral evolution of all transitions suggests 
that excluding systematic error for the PCA data to study 
the evolution of the disc normalisation in a relative sense is 
justified in practice. 

Occasionally, an observation that was adequately fitted 
when a 0.6% systematic error was added to the PCA data 
will converge on parameter limits or grossly unphysical val- 
ues when the data without systematic error are fitted. In 
these situations, the data without systematic error added 
are fitted by hand in order to achieve a fit with more realis- 
tic best-fit parameters. All errors presented in this work are 
la and are appropriately propagated when making relative 
measurements^ 

4-3.2 Disc Spectral Evolution: Variable / co i 



A variable colour correction 
from both a numer ical (e.g., 



iMerloni et al 



factor has been considered 
Shimura & Takaharal 1 19951 ; 



.. 2000 1; iDavis et all 2005) and observat onal 
(e.g.. |Punn et al.ll201ll ; iRevnolds fc Miller] 120111 ) approach, 
with 1.4 < / co i < 3 generally believed to be a conserva- 
tively broad allowable range, implying at most a factor of 
~ 2 change in / cc j. However, this range in / co i does not 
consider regimes extending to the low disc luminosities as- 
sociated with the HIS and LHS nor the effects of disc inho- 
mogeneities and strong magnetisation. Spectra that are less 
disc-dominated are harder; thus, requiring an increased / co i. 
The possibility of extending the upper / co i bound for discs 

6 Two identical data sets, one with a small systematic error added 
(X) and one without (Y), are nearly perfectly correlated (i.e., the 
correlation coefficient, px.Y = cov(X, V*)/(ffxffy) ~ 1). Since 
both X and Y only differ in the uncertainties associated with the 
data, their covariance is given by, cov(X , Y) = cov(X,X) = a x . 
The resulting correlation coefficient becomes, px.Y = ffx/ffy — 
1, assuming that the systematic errors added to Y are small, 
which is true for our purposes. 



© 0000 RAS, MNRAS 000, 000-000 





1 1 1 1 1 1 1 1 1 1 1 

; m= 0.0016 


i | i i i i 


5 1 




b — OB- Dh O - 





' , , , , 1 , , , , 1 , , 


, 1 , , , 



52380 
2P — r 



52390 



52400 



52410 



1 1 1 I 1 1 1 1 
m = -0.0043 

~OT> 



52700 



52720 



52740 



52760 



I 1 1 1 1 I 1 1 
m= 0.0019 



O- o- 



I i i_ 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 
- 00 D-O-OCP ®> - 



53190 53200 53210 53220 53230 53240 
2 



1 I 1 1 
m = 0.0021 



_L 



_L 



_L 



_L 



54120 



54130 



54140 



54150 



i , , i— | , , , , , , , , I 1 1 
m = -0.0035 



w- 



J_ 



54220 
2 



54230 



54240 



54250 



I 1 1 1 1 I 1 1 1 1 I 
m = -0.0007 

-Ate- 



_L 



J_ 



_L 



_L 



J_ 



55280 55290 55300 55310 55320 55330 
2 



* 1 
c 1 



'I 1 1 1 1 I 
m = -0.0016 



1 1 I 1 1 1 1 I 11 1 
* # 6 * ,«r- 



J I I I I l_ 



55570 



55580 55590 
MJD 



55600 



Spectral hardening in state transitions 
pow / bknpow comptt simpl 



11 



10 

1 

0.1 

10 

1 

0.1 

10 

1 

0.1 

10 

1 

0.1 

10 

1 

0.1 

10 

1 

0.1 

10 

1 

0.1 



I ' I . I . I ' I ' I ' + I . I ' I ' I ' I . I . + I ' I ' I ' I . I ' I 




I . I ■ I ■ I I I I I I I I ■ I I I I I I I ■ I ■ I I I I I I I I ■ I . I 



20 40 20 40 20 40 
t - 1 (Days) t - 1„ (Days) t - 1 (Days) 



Figure 6. Relative comparison of diskbb normalisations, Kbb rel 
(see text for definition), for each transition. The slope, m, from 
a linear fit to the data is displayed in the upper left corner for 
each plot. A zero-slope means that relative changes in A"bb with 
time are unaffected by neglecting to add systematic errors to the 
PCA data, indicating that the systematics remain constant from 
observation to observation for the time intervals considered. The 
flatness of -Kbb rel justifies the use of relative measurements for 
studying the disc evolution throughout a state transition. From 
top to bottom, the panels show the time evolution of -Kbb rel f° r 
transitions R02, D03, R04, R07, D07, R10, and Dll. Glyphs with- 
out error bars are larger than their respective la error. 



Figure 7. Evolution of the colour correction factor over all state 
transitions derived from best-fit diskpn normalisations where the 
inner disc radius was held at 6 i? g . The left, middle, and right 
panels correspond to fits using the pow/bknpow, comptt, and simpl 
continuum models, respectively. The abscissa shows time elapsed 
since the start time, to, of the transition (see Table [TJ. Error 
bars account for the maximum conceivable ranges of the mass, 
distance, and inclination of GX 339—4 (see text). The dashed lines 
enclose 1.4 < / co i < 3, which is the adopted allowable range in 
/ co l. The transition ID corresponding to each row is indicated 
on the rightmost ordinate. Within the large uncertainties, a non- 
truncated disc with a variable / col can adequately explain the 
observed disc spectral evolution. 
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in the HIS and LHS (i.e., non-disc-dominated regimes where 
idisc < 0.01 I/Edd) and whether a colour-corrected black- 
body is an appropriate description of the disc spectrum in 
this regime remain open questions. 

Using the diskpn model with the inner disc radius fixed 
ensures that the evolution of the disc normalisation is due 
entirely to a variable colour correction factor, which can 
be computed direc tl y from Equation [6] as demonstrated by 
iRevnolds fc Miller] ( 201lh . We elect to explore the degree 
of change in / co i alone that is required to explain the disc 
spectral evolution for a fixed inner disc radius of 6 R s . This 
choice for _Rj n is motivated by considering a disc extending 
down to the ISCO for a Schwarzschild black hole and is the 
minimum value allowed by diskpn. 

To measure the evolution of absolute / co i values required 
to explain the disc evolution in each transition for a fixed 
inner disc radius, we re-fit all of the spectra with 0.6% sys- 
tematic error added to the PCA data using the best-fit mod- 
els, but replacing the diskbb component with diskpn and 
freezing R[ n at 6 R s . Figure [7] shows the evolution in / co i 
that is necessary to account for each transition, with the er- 
ror bars incorporating the best-fit K pn error and the range 
in intrinsic parameters of GX 339-4. The values adopted 
for the GX 339-4 observables, along with their (lower, up- 
per) limits, are M = 7 (5.8, 12) M Q , D = 8 (6, 10) kpc, and 
i = 29° (10°, 60°) (see g2T]and references therein). Figure [7] 
shows that, within the uncertainties, the disc spectral evolu- 
tion can potentially be explained by a variable 1.4 < / co i < 3 
for an inner disc radius that remains at 6 R s throughout each 
state transition. 

Studying the evolution of absolute f co \ values for a fixed 
Rin is limited by the large uncertainties in the mass, dis- 
tance, and inclination of GX 339-4. Given that an inner 
disc that remains at 6 R s combined with a modestly variable 
/coi may be a viable model for accretion disc evolution in the 
state transitions studied here, we use the relative measure- 
ment technique introduced in £14.3. 1 1 to better understand 
the degree of change in / cc j that is required to explain the 
observations without appealing to disc truncation. 

The PCA data are replaced with a duplicate copy where 
no systematic errors have been added. Next, the diskpn fits 
where R[ n = 6 R s are re-fitted and a best-fit is found for 
each observation. This is done for each of the three Comp- 
tonisation models pow/bknpow, simpl, and comptt. For each 
observation in a given transition, / co i is computed from K pn 
using Equation [6] without specifying the intrinsic parame- 
ters (i.e., M, D, i). The resulting / co i values are normalised 
to the / co i from a representative HSS observation, / co i,hss, 
and furt her normalised to the adopted / co i = 1.4 value for 
the HSS l|Davis et al.ll2005h . yielding the relative colour cor- 
rection factor, 



fc, 



1.4 



/col 



fc 



(8) 



Supposing that the inner disc remains fixed at 6 R s 
throughout each state transition, Figures IBl llOl show the rel- 
ative change in / co i alone required to explain the observed 
disc spectral evolution. Notably, the pow/bknpow model re- 
quires more drastic change in / co i compared to comptt and 
simpl. Under the presumption that / co l,rel = 3 is a conser- 
vative upper bound, the notion of a variable / co i alone is 
perhaps possible for transitions R02, R07, R10, and Dll; 
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Figure 8. Relative measurements of the evolution of / co i for each 
transition derived from fits using pow/bknpow to model Compton- 
isation. The two lower dashed lines enclose 1.4 < / co i < 3, while 
the upper dashed line marks / co i = 5. The / co i re i measurements 
for each transition were made relative to the filled glyphs and 
normalised to / co i = 1.4, which is appropriate for the HSS. The 
disc evolution can be explained for transitions R02, R07, R10, 
and Dll by a non-truncated disc with 1.4 < / co i,rcl < 5. 



however, an evolving / co i is an inadequate description of the 
evolving disc spectrum for transitions D03, R04, and D07, 
regardless of the Comptonisation model employed. For the 
later set of transitions, one must invoke disc truncation or an 
alternative disc model to satisfacto r ily ex plain the observed 
disc evolution. IRevnolds fc Milled { 201 lh speculate that if 
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Figure 9. Same as Figure [8] but with the Comptonisation model 
comptt replacing pow/bknpow. The required evolution in / GO i re i is 
systematically reduced compared to results from the pow/bknpow 
continuum model. 



Figure 10. Same as Figure \E\ but with the Comptonisation 
model simpl replacing pow/bknpow. The required evolution in 
/col rel is systematically reduced compared to results from the 
pow/bknpow continuum model. 



the colour-corrected blackbody is a viable model of the LHS, 
then /coi may be as large as ~ 5 based on spectral fits to a 
comprehensive set of BHB spectra. Inspecting Figures [91101 
where the comptt and simpl models were adopted, all seven 
transitions essentially conform to this extended / 0O i range; 
although, making definitive claims regarding the physically 
plausible / co i range and applicability of the MCD model in 
the LHS is beyond our scope. 



4.4 A Hot Inner Disc Can Masquerade as an 
Increased / co i 

Observational evidence suggests that GX 339-4 harbors 
a near-maximally spinning black hole (see t|2.ip . Conse- 
quently, the inner disc extends deep into the black hole 
potential and relativistic effects become non-negligible, ren- 
dering Newtonian disc models inappropriate. The disc mod- 
els used in this work to derive / co i measurements impose 
an effective lower limit on the inner disc radius of 6 R K , 
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whereas the ISCO of a Kerr black hole with a, = 0.998 
lies at 1.24 _R g , corresponding to a substantially hotter in- 
ner disc. Consider attempting to use a Newtonian prescrip- 
tion to model an accretion disc around a Kerr black hole, 
where R[ n = f?isco- The Newtonian disc model is restricted 
by i?i n > 6 R g and will attempt to fit the true hotter 
disc, located at -Risco < 6 R s , by increasing the colour 
correction factor as follows. The colour temperature deter- 
mined by the model depends on both f co \ and T c g accord- 
ing to T C oi = /coiTeft- The maximum effective temperature 
of the Kerr disc exceeds that of the Newtonian model with 
i?in = 6i? g ; therefore, the only way for the model to achieve 
the colour temperature demanded by the data is to increase 
/coi to an erroneously high value. 

Using the XSPEC command fake it none, we generate 
a suite of simulated disc spectra with the kerrbb model of 
a spinning black hole with an RXTE response matrix and 
integration time of 2000 seconds. The cases considered are a 
Schwarzschild black hole (a* = 0, f] = 0.057) and m aximal 
Kerr black hole (a, = 0.998, r) = 0.30) (|Thornelll97i ). each 
having mass 7 Mq located at a distance of 8 kpc with a 
disc inclination of 30°. The inner disc torque is set to zero 
and the effects of both self-irradiation and limb-darkening 
on the disc spectrum are switched on. We pair mass accre- 
tion rates with colour correction factors that are reasonably 
representative of the spectral state corresponding to each M 
that we select. For the a, = black hole, we simulate disc 
spectra for M = 0.05 MEdd with / co i = 1.4, attempting to 
mimic a HSS/SIS accretion disc; and M = 0.005 MEdd with 
/coi = 3, which is more representative of a HIS/LHS obser- 
vation. Similarly, for the a* = 0.998 black hole, we choose 
M = 0.01 M E dd with /coi = 1.4 and M = 0.001 M E dd 
with /coi = 3. Extending M beyond these ranges for a fixed 
/coi resulted in either simulated data with unconstrained 
uncertainties or unacceptable fits (i.e., xl 2> 1) using the 
procedure described below. 

Each spectrum is fit over the energy range 2.8-10 keV 
using the diskpn model with Ri n fixed at 6 R s . The best- fit 
colour correction factor is computed from the normalisa- 
tion using Equation [6] Table [2] presents the results of this 
exercise. In general, the psuedo-Newtonian model, diskpn, 
reliably recovers the true / co i value for a non-spinning black 
hole, but requires a factor of ~ 1.9 increase in / co i for a 
maximally spinning black hole. This demonstrates that the 
inability of Newtonian disc models to properly model Ri n 
and T a g for Kerr black holes can result in unphysically large 
/coi measurements to account for the hotter inner disc. 

One may wonder if relative measurements of / co i de- 
rived from Newtonian disc models, as done in this work, are 
affected by the inability of these models to properly repre- 
sent the hotter discs around spinning black holes. A colour 
correction factor evolving in a HSS —> LHS transition from 
/coi = 1.4 — ¥ 3 corresponds to a factor of 2.1 increase in 
/coi- Inspection of Table [2] indicates that the diskpn model 
requires a factor of ~ 2.1 increase in / co i for an order of mag- 
nitude decrease in M for both the Schwarzschild and max- 
imal Kerr black holes. This result is fully consistent with 
the 'true' / co i evolution and shows that relative measure- 
ments of /coi made with Newtonian disc models are reliable, 
even when the black hole under consideration is maximally 
spinning. 



a* 


M 


^col,max 




J col, true 


/col 


~Kv 




(M Edd ) 


(kcV) 













0.05 


0.416(2) 


0.177(8) 


1.4 


1.39(2) 


1.35 





0.005 


0.493(5) 


0.0098(8) 


3.0 


2.87(6) 


1.26 


0.998 


0.01 


0.455(5) 


0.015(2) 


1.4 


2.58(9) 


0.87 


0.998 


0.001 


0.56(1) 


0.0006(1) 


3.0 


5.8(2) 


0.71 



Table 2. Results from fitting simulated spectra generated from 
the kerrbb model of a Kerr black hole with the pseudo-Newtonian 
model, diskpn, where the inner disc radius was set to 6 R g . The 
first two columns are the black hole spin and mass accretion rate 
in Eddington units used to generate the simulated disc spectra. 
Proceeding from left to right, the columns are the best-fit maxi- 
mum disc colour temperature, best-fit disc normalisation, colour 
correction factor used to generate the disc spectrum, colour cor- 
rection factor recovered from the best-fit normalisation, and re- 
duced x 2 . All models had 15 degrees of freedom. Uncertainties on 
the last significant digit are given in parentheses and correspond 
to the 90% level. 



4.5 Disc Temperature and Luminosity 

The luminosity of a geometrically thin, optically thick disc 
is theoretically expected to obey a relation of the form 
idisc = ffATeff, where A is the radiating surface area. The 
disc luminosity depends on both A, which scales as Rf n , and 
effective temperature, which depends on / co i. The purpose of 
this section is to explore the Ldisc — RfnT^ff relation, partic- 
ularly extended to low luminosities (i.e., the HIS and LHS), 
and the effect of introducing the possibility of a variable 
colour correction factor. 

For the ezdiskbb disc model of a N ewtonian disc with 
a zero -torque inner boundary condition, I Zimmerman et all 
(2005) finds the disc luminosity is given by, 



Ldisc = 73.9<t 



( Tcol 



V u 



Rin 



(9) 



Observationally, we measure an unabsorbed disc flux, 

In practice, Fdi sc is obtained by operating the XSPEC cf lux 
command on the disc component with the normalisation 
frozen to its best-fit value. Unabsorbed disc fluxes and unab- 
sorbed total model fluxes, F to tai, are calculated over 0.1-10 
keV and 0.1-100 keV, respectively. Figure [TT1 shows the re- 
lation between unabsorbed disc flux and maximum effective 
temperature for each of the three Comptonisation models 
pow/bknpow, comptt, and simpl. T co i,max is measured from 
fits to all observations with the disc component modelled by 
ezdiskbb and T c ff,max = Tcoi,max//coi, assuming a constant 
/coi = 1.4. The disc fractions, defined by Equation [7] are 
given by the glyph colours. We note that this work includes 
observations with very low disc fractions compared to pre- 
vious investigations of the luminosity-temperature relation 
for BHBs, whi ch tend to be restricted to the HSS where 
DF > 0.8 (e.g., lGierliriski fc Dondl2004l ). 

If one assumes a constant disk geometry, the disk emit- 
ting area is unchanging and we expect to recover the scaling, 
-fdisc oc T e 4 ff , denoted by the diagonal lines in Figure [TT1 for 
various choices of Ri n and / co i- While the theoretical rela- 
tion is achieved for disc luminosities, Ldisc > 0.05I/Edd, clear 
departure is evident below this threshold where the BHB is 
categorized in the HIS or LHS. 
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Figure 11. Unabsorbed disc flux (0.1-10 keV) vs. maximum disc 
effective temperature (T cffjmax = T col>max // col ) for all transi- 
tions, where T co i max is measured using ezdiskbb to model the 
disc component and a constant / co i = 1.4 is assumed. Comp- 
tonisation is modelled using the XSPEC models pow/bknpow (top 
panel), slmpl (center panel), and comptt (bottom panel). From 
top to bottom, the horizontal dashed lines indicate idisc/^Edd = 
10 _1 , 10~ 2 , 10~ 3 , where L 

disc was converted to Fdi sc using Equa- 
tion [10] with D = 8 kpc. The Eddington luminosity is given by, 
L B dd = 4rrGMmpc/ffT, where m p is the proton mass, ax is the 
Thomson cross section, and we adopt M = 7 Mq. Uncertain- 
ties in M and D are not incorporated because the shape of the 
^disc - T e g max relation is insensitive to these intrinsic parameters. 
The diagonal dashed lines represent the expected Fdi sc oc max 
relation for an optically thick, geometrically thin, zero-torque disc 
(Equation |9| assuming the aforementioned intrinsic parameters 
and constant i?; n = 1 R e (black line), 6 R g (red line), 20 R g 
(blue line), and 100 R g (green line). Glyph colours denote the 
disc fraction of each observation. The key in the bottom right 
identifies the glyph associated with each transition ID. When a 
constant / co i is assumed, departure from the Ldisc - T^S ' s ev ident 
for low-luminosity discs characteristic of the HIS and LHS. 



The diskbb and ezdiskbb models are purely Newto- 
nian and correction factors accounting for various general 
relativistic and vi ewing; angle effects (|Cunninghaml Il975l ; 
IZhang et all 1 19971 ) have been applied in previous obser- 
vational stu dies of the luminosity-temperature relatio n for 
BHBs (e.g., Gierliriski fc Dondl200i iDunn et al.ll20Tll ). Se- 



Figure 12. Relative measurements of the colour correction factor 
shown against the maximum colour temperature of the disc from 
best-fits using the XSPEC Comptonisation models pow/bknpow (top 
panel), simpl (center panel), and comptt (bottom panel). Results 
from all transitions are shown together, with the disc fraction 
given by the glyph colour. The key in the upper left portion of 
the top panel identifies the glyph associated with each transition 
ID. T ro 

l,max is obtained from the best-fit ezdiskbb parameter 
and / co i,rcl is computed using the same procedure as in Figure [8] 
(see text). Horizontal dashed lines mark / CD i Ie l = 1.4,3,5. For a 
disc inner radius fixed at 6-R g , the majority of observations yield 
1-4 < /col rel 5 an d show a roughly linear relation between / co i 
and T col . 



lecting the proper correction factors requires knowledge of 
the disc inclination and spin, for which there is considerable 
debate for GX 339-4 (see £|2. 1[) . Fortunately, the correction 
factors are independent of flux and temperature; therefore, 
we can safely neglect incorporating these modifications with- 
out altering the shape of the -Fdi sc — Toff diagram, irrespective 
of the accuracy of the Fdisc and T e g values. 

The effective temperature plotted in Figure [TT] is calcu- 
lated under the assumption of a constant colour correction 
factor, which will subsequently be relaxed. Deriving relative 
measurements for the colour correction factor according to 
Equation [H] and as described in H4.3.21 Figure [T2"l shows the 
relationship between / co i and the maximum colour temper- 
ature, which is a parameter in the disc models. Supposing 
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(keV) 



Figure 13. Unabsorbed disc flux (0.1-10 keV) vs. maximum disc 
effective temperature for all transitions, where T co i max is mea- 
sured using ezdiskbb. Relaxing the assumption of a constant / co i, 
all observations were re-fitted with ezdiskbb replaced by diskpn 
with _Ri n = 6 R g . Relative values of / co i were computed from 
the best-fit diskpn normalisation, using the same procedure as in 
FigureH The abscissa shows T offjmax = T col , max // col , rcl . Glyph 
colours correspond to the / co l,rcl of each observation. The various 
lines are the same as those in Figure [TTJ The L<ji sc <x relation 
is recovered for disc luminosities, L disc > 10 — 3 I/Eddi when / col 
is considered to be variable. 



that the inner disc remains fixed at R ln = 6 i? g , which is 
the implicit assumption behind all of our / co i,rei measure- 
ments, the colour correction factor increases with colour 
temperature approximately linearly. The majority of obser- 
vations yield colour corrections that fal l within the range 
1.4 < /col.rel < 5, in agreement with iRevnolds fc Miller] 
(|201lh . Notably, the models simpl and comptt, which are 
rooted in physical considerations of Comptonisation, con- 
strict the range of / co i,rci as compared to the phenomeno- 
logical and widely used powerlaw and bknpower models. Par- 
ticularly for best-fit colour temperatures, T co i imax > 1 keV, 
the powerlaw and bknpower continuum models may return 
erroneously low values for the disc normalisation. 

We attempt to explore whether the luminosity- 
temperature relation can be extended to lower luminosities 
characteristic of the HIS and LHS by allowing for a variable 




40 

t- 1 (Days) 



Figure 14. Time evolution, as measured from the start of the 
transition listed in Table[JJ of the inner disc radius computed from 
the ezdiskbb and diskpn normalisations (Equation [TT} when Ri n 
was fixed in the diskpn model at 6 R g . The dashed line marks 
Bi n = 6 -Rg . The key in the upper right corner denotes the glyphs 
corresponding to each transition. The constancy of i?; n implies 
that relative measurements made from combining results from 
fits made with ezdiskbb and diskpn arc unaffected by the minor 
differences between those disc models. 



/coi- Figure [T3l shows the -Fdi sc — ^eg diagram with the only 
difference to Figure [TT] being that the effective temperature 
is now calculated by, T cf f, max = T co i, m ax//coi,roi. In other 
words, Figure [TTJ is recomputed using the / co i,rci values of 
Figure [T2l in place of the constant / co i = 1.4. Accepting the 
dynamic range in f co \ suggested by Figure 1121 Jdisc oc T^ s 
is recovered for disc luminosities, -Ldisc > 10 -3 I/Edd- Fur- 
thermore, the relation falls neatly on the expected track cor- 
responding to the emitting surface area remaining constant 
with the inner disc located at 7?i n = 6 R s . 

Caution must be exercised when combining different 
disc models. The relative measurements / co i,rei used to re- 
compute TcflSmax in Figure [T3] called upon diskpn and 
ezdiskbb separately to model the disc component in order 
to take advantage of the strengths of each model for our 
purposes. One may reasonably ask if these two disc models, 
which both adopt the zero-torque inner boundary condition, 
are consistent. The significant difference between the models 
is that diskpn approximates the general relativistic effects 
on the accretion disc with a pseudo-Newtonian potential, 
while ezdiskbb is strictly Newtonian; therefore, discrepan- 
cies between the two may materialize in the spectral fits. 
Combining the normalisations from each disc model, given 
by Equations [S] and [5J the inner disc radius is obtained, 

1/2 



Rin 

Ra 



= 0.0677 



(11) 



Since Ri n is a separate parameter in diskpn instead of being 
incorporated into the normalisation, we expect to exactly 
recover this value for i?i n from Equation [TTJ if diskpn and 
ezdiskbb are identical. Figure [14] shows R ln computed from 
Equation [TT] for fits to all observations where R ln = 6 i? g in 
the diskpn fits. Values of R ln falling on 6 i? g would indicate 
that the inner radius enforced in the diskpn fits was exactly 
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recovered in the independent ezdiskbb fits, while deviations 
from 6 R g arise from differences in the ezdiskbb and diskpn 
models. Although the inner disc radius value of 6 R g set 
in diskpn is not exactly recovered, Ri n remains constant 
throughout each transition at ~ 7.5 R s , merely a systematic 
offset that does not affect the interpretation of Figure If 31 

4.6 Low Frequency Quasi-Periodic Oscillations 

LFQPOs are observed in the temporal power spectra of 
BHBs and remain a long-standing puzzle of active ob- 
servational and theoretical interest. Explanations for the 
origin of LFQPOs appeal to an exhausti ve array o f pro- 
duction models inc luding trapped wav es (lKatdll990l ). dis- 
coseismic modes dNowak fc Wagoner] Il99ll . Il992l . Il993l; 
Reynolds fc Miller! \200$l Lense-Thirrin g prec ession (|lpserl 
19961: IStella et al.1 Il999l; llngram et al.l I2009T). accretion- 



eject i on instability (I Tagger fc Pellatl 19991 ; IVarniere et al.l 



I2OO2I ; IVarniere fc Taggerj^OOd) , an d a truncated disc LHS 

it 2 



geometry (jGiannios fc Spruitll2004l ). Recent simulations of 
global, magnetised accretion discs suggest that LFQPOs 
have an origin in th e corona rather than the disc directly 
(|Q'Neill et al.ll201ll ). While the origin of LFQPOs is not 
firmly established, the details of the disc geometry may pro- 
vide clues to their nature. 

X-ray spectra alone cannot distinguish between a trun- 
cated disc and changing disc structure/coronal activity as 
the underlying accretion geometry. Given a plausible pro- 
duction mechanism, LFQPOs have the potential to be di- 
agnostics of the inner acc retion flow. In the LFQPO excita- 
tion model envisioned bv lGiannios fc Spruij |2004 ). quasi- 
spherical oscillations of the RIAF interior to a truncated 
disc in the LHS produce modes with frequencies near the 
Keplerian frequency of the inner disc, 

-2/3 

(12) 



Rkc 

Hey 



(GM 

y— 2 ™Qpo 



where z^qpq is the frequency of the LFQPO. In this frame- 
work, iGiannios fc Spruid (|2004f ) require that the RIAF, and 
hence the recessed inner disc, extend out to Ri n ~ 100 R g . 

LFQPOs observed in GX 339-4 and their frequency 
migration over the course of the state transitions studied 
here has received considerable attention in the literature. 
The frequency of LFQPOs is observed to decrease (increase) 
over a decay (rise) transition with a typical range 0.1 Hz < 
vqpo < 10 Hz, implying that accretion rate and QPO fre- 
quency are intricately linked. In an attempt to test the trun- 
cated disc LFQPO production model, which predicts that 
J?Ke P — Rin > 100 i?g in the LHS, we collect measured 
vqpo values for all of the GX 339-4 state transitions consid- 
ered dBelloni et al.l 120051; IShaposhnikov fc Titarchuk) |2009| ; 



IStiele et al.l l201ll ; iMotta et al.l l201ll ; iNandi et al.l I2OI2T ). 
Given the uncertainty in black hole mass for GX 339-4, we 
seek relative measurements of -Rk cp normalized to the ISCO 
for a non-spinning black hole in the HSS, 



R 



Kcp.rcl 



6R S 



Rkc 



(13) 



^Kcp.HSS, 

where -Rk cp ,hss is the Keplerian radius computed for a given 
transition from the j^qpo of the observation nearest to the 
HSS in time using Equation ll2l without specifying M. Sup- 
posing that i?Kc P is associated with R[ n , Figure [15] shows 
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Figure 15. Relative measurements of the evolution of -RKep,rel 
for each transition (see text for definition) derived from LFQPO 
frequency evolution. The dashed lines mark R ln = 6,20,100 R g . 
The i?Kep,rel measurements for each transition were made rela- 
tive to the filled glyphs and normalized to Hi n = 6 -R g , which is 
appropriate for a Schwarzschild black hole in the HSS. 



the evolution in disc inner radius implied by the observed 
migrating LFQPOs for each transition. With the exception 
of transition Dll, the Keplerian frequencies implied by the 
LFQPOs reach out to ~ 50 - 100 R g in the LHS, which 
is broadly consistent with the degree of recession required 
by the RIAF/truncated disc LFQPO excitation model. If 
LFQPOs are associated with the inner disc in a simi lar man- 
ner to that described bv lGiannios fc Spruij l |2004 ). then a 
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truncated disc geometry may be a plausible description for 
state transitions. However, the spectral fits in this work sug- 
gest that Rin does not truncate to such extreme disc loca- 
tions, which would demand reconciling. 



5 DISCUSSION 

We have shown that a moderately increasing f co \ as the disc 
fraction decreases can adequately and consistently describe 
the spectral evolution of the accretion disc in GX 339-4 over 
multiple state transitions. The purpose of this work is not 
to debunk a truncated disc geometry as a viable model for 
the LHS of BHBs. We recognize that the observations in 
this work interpreted with simplistic disc models are unable 
to distinguish between changes in Ri n and/or / co i. Rather, 
we demonstrate that reasonable variations in / co i provide an 
alternative explanation to disc truncation for accretion disc 
evolution. Motivated by this result and supposing that f co \ 
is variable, we consider the impacts that this would have on 
the interpretation of BHB observations, ultraluminous X- 
ray sources, and models of the accretion flow across spectral 
states. We then discuss the physically expected range of at- 
tainable values for / co i and conclude the discussion with a 
cautionary commentary on the widely used disc models. 



5.1 Implications of a Variable Colour Correction 

Imposing an unchanging inner accretion disc radius located 
at the ISCO of a Schwarzschild black hole, we measured the 
relative changes in / co j necessary to adequately fit RXTE 
spectra of GX 339-4 over seven state transitions. Arbitrarily 
selecting a normalisation of / co i = 1.4 for the HSS and de- 
noting measurements of / co i relative to this value by / co i,rei 
(see Equation [8]), we tracked the evolution of / co i,rel for each 
transition (see Figures [8TT0)> . Transitions R02, R07, RIO, 
and Dll require colour corrections conforming to the range 
1-4 ^ /coi.rci i$ 3, while transition R04 demands the range 
be extended to 1.4 < / co l,rel i$ 5. These ranges in / co i evolu- 
tion, which can account for the majority of the transitions, 
correspond to a factor of ~ 2.0 — 3.5 increase in / co i for a 
HSS — > LHS transition. Transitions D03 and D07 require 
more extreme evolution in / co i,rci, with values as high as 
/coi,rei ~ 10. Our main finding is to demonstrate that an 
evolving / co i provides an alternative and adequate descrip- 
tion of the changing disc spectrum during a state transition, 
which is perhaps a more attractive scenario to imagine than 
the sudden evacuation or replenishment of the inner accre- 
tion flow. 



5.1.1 Recovering L — T e g at Low Luminosities 

A blackbody source of constant emitting area is expected to 
obey the L — T* ff scaling, which relates the source luminosity 
to its effective temperature. Accretion flows corresponding 
to the HSS of BHBs, where the thermal accretion disc dom- 
inates the total flux and the inner disc location does not 
devia te from the ISCO, nicely follow this expected relation 
(e.g.. iGierlinski fe~D onc 2004) ; however , deviations arise in 
the HIS and LHS fe.g jDunn et al.ll201ll ). This work encom- 



finding that a f co \ that evolves with spectral state can re- 
produce the theoretical Ldisc — T^g across all spectral states. 
Studying the disc-dominated states (i.e., DF > 0.8) of many 
BH Bs over a typical luminosit y range, 0.01 < -Ldisc/^Edd < 
0.5, IGierlinski fc Pond |2004l ) found that the L d isc - T* H 
relation holds for a constant f co \. We observe similar and 
consistent behaviour for GX 339-4 in the HSS. As the disc 
fraction drops and the source enters the HIS or LHS, signif- 
icant departures fr om Ldi BC — T^ are apparent in Figure [TTJ 



passes the disc luminosity range, KT 3 < Ldisc/iEdd < 0.5, 



in agreement with lDunn et all ( 201ll ). These 'spurs' on the 
luminosity-temperature diagram can be interpreted as either 
a decrease in the inner disc radius, and hence the emitting 
area, or the result of neglecting to account for a decrease 
in T e s = T C oi//coi caused by an increased / co i in the LHS. 
Of cour se, some combinat ion of these two effects is possible 
as well. iDunn et all (|201lh measured the incremental degree 
of change in / co i from a canonical HSS value that was re- 
quired to return the spurs to the expected I/disc <x T^ g law 
for a constant emitting area (see their Figure A.l). For GX 
339-4, they found that fairly modest changes in / co i, consis- 
tent with the relative changes we determined, could restore 
the theoretical relation; however, they argue that the most 
severe departures from Z/disc — T^ s , which occur for the low- 
est disc fractions, are not likely the result of a variable f co \. 
Here, we suggest that the physically realizable range of / co i 
is perhaps not as narrow as commonly assumed (see i]5.2[) 
and that a thin disc geometry with a fixed inner radius and 
evolving / co j is a viable phenomenological model for BHB 

state transitions. 

We expand on the work of lDunn et al.l (|201ll ) by break- 
ing the observations of a BHB into individual transitions, ex- 
ploiting relative measurements, and exploring empirical and 
physical Comptonisation models. Comparing Figures 1 81101 
we find that employing physically motivated Comptonisa- 
tion models that avoid the 'flux stealing' phenomenon from 
the disc component require less severe / co i evolution for any 
given transition. Inspection of the panels of Figure [13] shows 
that the Ldi sc — T* s relation is recovered when a changing / co i 
is incorporated, regardless of the choice of Comptonisation 
model. Relative measurements of / co i broken down by tran- 
sition in Figures [81 101 hint at an underlying difference in the 
details of the accretion flow depending on whether the source 
is in the rise or decay stage of outburst. Qualitatively, the 
decay transitions demand more extreme and perhaps phys- 
ically unattainable / co i evolution, while the rise transitions 
fit more comfortably into the idea of a moderately variable 
/coi. Given that decoupling the true behaviour of the inner 
disc from that of the colour correction factor is impossible 
with simplistic disc models and RXTE data, we leave the 
possibility of different accretion scenarios between rise and 
decay transitions as pure speculation. 



5.1.2 The Ultraluminous State 

Adopting a variable colour correction factor parameterises 
the various physical processes in BHB accretion into a single 
factor. This work suggests that moderate changes in / co i may 
conceivably account for the wide range of spectral states ob- 
served when a BHB undergoes a state transition. Recently, a 
new spectral class, dubbed the ultraluminous state (ULS), 
has been postul ated for the class of u ltraluminous X-ray 
sources (ULXs) (|Gladstone et al.l l2009h . Here, we explore 
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whether the ULS, which has never been observed in a con- 
firmed BHB, is consistent with the interpretation of spectral 
states arising from physical processes that are observation- 
ally manifested as an evolving colour correction factor. 

ULXs are X-ray sources that exceed the Eddington lu- 



^disc ' ^total 



minosity of stellar mass black holes (Lx > 10 erg s 



1 ), but 



are not spatially coincident with galactic nuclei; therefore, 
ULXs are distinct from AGNs. Two possible explanations 
for ULXs are: (1) intermediate mass black holes (IMBHs; 
M ~ 10 2 — 10 4 Mq) accreting at sub- Eddington rates, but 
the production rate of binaries composed of an IMBH and 
donor star is li kely too low to account fo r the number of 
observed ULXs |Madhusudhan et alj|2006h : and/or (2) stel- 
lar mass black holes accreting at super-Eddington rates, 
which requires a mechanism for achieving super-Eddington 
accretion, such as beaming /viewing effec ts of outflows (e.g., 
iKing et al.l l200ll; iBegelman et alTl2006l ) or disc inhomo- 
geneities (lBegelmanll200ll , |2002| . |2006h . 

The spectral char acteristics of ULXs were studied by 
IGladstone et al.l (|2009l ), leading them to propose that ULXs 
are super-Eddington accreting stellar mass black holes with 
ubiquitous spectral features. These commonalities suggest 
that ULXs are indeed a spectral class of BHBs, dubbed 
the ultraluminous state (ULS); however, this accretion state 
does not appear to arise in the well-studied Galactic BHBs. 
The ULS is characterised by both a broadened, hardened 
thermal disc component and a break in the power-law com- 
ponent at lower energies (> 3 keV) than typically observed 
in less luminous BHB states. The IGladstone et al.l (|2009T ) 
ULX sample was adequately fitted with a physically moti- 
vated disc plus Comptonising corona model, where the best- 
fit models called for optically thick coronae (r p ~ 5 — 30) and 
cool disc colour temperatures (T co i jma x ~ 0.2 keV). Physi- 
cally, the broadened/hardened disc and low-energy spectral 
break may be attributed to a complicated inner geometry 
where an optically thick corona obscures and alters the en- 
ergetics of the inner disc. Normally, disc temperatures this 
low would imply a black hole of mass ~ 10 3 M Q (e.g., 
iKaaret et al . 2003); however, the measured disc temperature 
may be lowered if a significant fracti on of the accretion en- 
ergy is dissipated in the corona (e.g.. [Svensson &: Zdziarskl 
ll994T ). Based on the requirement of an optically thick corona, 
as opposed to the r p ~ 1 coronae typically observed in 
BHBs, sub-Eddingtion accretion onto an IMBH was ruled 
out. 

Figure [16] shows the coronal optical depth versus the 
maximum disc colour temperature for GX 339-4 deter- 
mined from the best-fits to the disc plus Comptonising 
corona XSPEC model, phabsx (comptt+diskpn) , with 7?i n = 
6 R s . Replacing the disc component with either diskbb 
or ezdiskbb produces essentially indisting uishable results. 
For co mparison, the results from Table 8 of lGladstone et all 
(2009) are shown for a sample of 12 ULXs determined from 
fits to high-quality XMM-Newton data, where diskpn and 
comptt were also adopted as the disc and Comptonising 
models, respectively!!] When GX 339-4 is in the HIS/LHS 



7 IGladstone et all J2009h do not specify the i?; n value used in 
the diskpn model. Our results suggest that the coronal properties 
and disc temperature are largely unchanged in the spectral fits, 
regardless of what one chooses for Ri n . 
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Figure 16. Coronal optical depth vs. maximum disc colour tem- 
perature for all transitions, where r p and T co i max are measured 
using comptt and diskpn with R ln = 6 R g , respectively. The disc 
fraction is given by t he glyph colour. The re sults from a sample of 
12 ULXs (Table 8 of lGladstone et alj 1120091 )1 are overlaid in black 
crosshairs. The HIS/LHS GX 339-4 RXTE data analysed in this 
work permit moderately optically thick coronae in the same range 
as found for ULXs. 



(i.e., low disc fractions), the typical best-fit coronal optical 
depth falls in the range, r p ~ 2 — 10. The coronal optical 
depth of the ULX sample is described by r p ~ 5 — 15 (r p ~ 



10 — 100) for discs with temperatures T cl 



< 0.5 keV 



(T c< 



0.5 keV). Figure [TBI challenges the claim that 



spectral fits to sub-Eddington BHBs do not permit moder- 
ately optically thick coronae. 

Although the data are consistent with an optically thick 
corona model, here we point out a strong parameter degen- 
eracy between the colour correction factor and the coronal 
optical depth. Of the three models we considered to repre- 
sent Comptonisation, all of which provided equally good fits 
to the data, only comptt incorporates the coronal tempera- 
ture and optical depth as physical parameters. Inspection of 
Figure[T2]reveals that the pow/bkn and simpl models require 
higher colour correction factors than comptt. Particularly for 
the HIS /LHS observations of GX 339-4, introducing a phys- 
ical corona drives / co i to lower values with preference for an 
optically thick corona (see Figure [T2")) . Therefore, the data 
cannot distinguish between an increased colour correction 
factor and an optically thick corona. 

In the IMBH model, ULXs are massive scaled up ver- 
sions of BHBs in the LHS (i.e., IMBHs accreting at sub- 
E ddington rates). The IM BH model for ULXs was discarded 
bv lGladstone et al.l (2009) because the measured coronal op- 
tical depths were larger than commonly observed in BHBs. 
In the case of GX 339-4 in the HIS/LHS, using the same 
disc plus Comptonisation model, we see that the best-fit 
coronal optical depths lie precisely in the same moderately 
optically thick range as the ULX sample of lGladstone et ahl 
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Figure 17. Coronal optical depth vs. coronal temperature, both 
measured using comptt. The disk model was diskpn with i?; n = 
6 R g . The disc fraction is given by the glyph colour. The result s 
from a sample of 12 ULXs (Table 8 of ^Gladstone et al.l ([2009)) 
are overlaid in black crosshairs. 



(2009). Therefore, the data do permit optically thick coronae 
in the HIS/LHS of GX 339-4, but we stress that the mea- 
surement of coronal optical depth from basic disc + corona 
modelling cannot be made definitive due to the degeneracy 
between r p and / co i. This possibility alleviates the hesitation 
for supposing that ULXs are analogous to scaled up BHBs 
on the grounds of coronal optical depth. Allowing for a vari- 
able colour correction factor, the broadened/hardened ULX 
disc spectra and low measured ULX disc temperatures nat- 
urally fit into the IMBH description; therefore, introducing 
arguments for reducing the disk temperature by dissipat- 
ing accretion energy into the corona become unnecessary. 
An important discrepancy in coronal properties between the 
HIS/LHS GX 339-4 data and the ULX sample is the tem- 
perature of the corona. While the data permit r p ~ 5 — 10 
for both the HIS/LHS and ULXs (see Figure Ho) , the asso- 
ciated coronal temperature for ULXs (T p ~ 1 — 3 keV) is an 
order of magnitude below that of GX 339-4 in the HIS/LHS 
(T p ~ 10 - 30 keV) (see Figure \T7\. We stress that the pur- 
pose of this section is to highlight an interesting possibility 
that the ULS may not be an intrinsic feature of BHBs; how- 
ever, this is highly speculative and further investigations are 
required to shed more light on the nature of ULXs. 

5.1.3 The Role of f co i in State Transitions 

A triggering mechanism for st a te tra nsitions is not yet well 
understood. ISchnittman et all (|2012f l demonstrated that all 
of the spectral states of BHBs could be recovered by vary- 
ing only the mass accretion rate, M, (i.e., effectively chang- 
ing the location of the disc photosphere, or the disc/corona 
boundary) in the magnetohy drodynamic genera l relativis- 
tic global disc simulation of iNoble et all l|2010l ). The disc 



was assumed to produce a diluted blackbody spectrum 
at the photosphere with a constant f co \ = 1.8, which 
could be further altered (e.g., broadened/hardened) by 
bremsstrahlung, synchrotron, and inverse Compton pro- 
cesses in the corona. The emergent spectrum was com- 
puted with a post-processing radiative transfer technique 
for choices of M representative of the observed BHB states. 
In this simulation, the inner accretion disc remained ex- 
tended down toward the ISCO. Only details of the coronal 
physics for a given location of the disc photosphere (i.e., for 
a given value of M) provided alterations to a diluted black- 
body spectrum and all BHB spectral states were reproduced 
within this framework. This result emphasises the impor- 
tance of coronal physics in modifying the underlying disc 
spectrum, suggesting that perhaps / co i is a variable func- 
tion of M and can be dependent on processes external to 
the disc. 

A changing mass accretion rate plays a key role in state 
transitions, but is likely not the whole picture, nor is it clear 
what mechanism drives changes in M. For an inner disc 
remaining at the ISCO during a state transition, a large 
change in M is implied. In the HSS, the inner portions of 
the disc are expected to be radiation pressure-dominated, 
where instabilities may le ad to an inhomo geneous disc struc- 
ture and a low / co i ~ 1 (Bcgelman 2006). As M decreases, 
the ratio of radiation-to-gas pressures drops and below some 
idisc/i/Edd the density inhomogeneities would disappear and 
/coi would increase. If the di sc is strongly magnetised, v alues 
°f /coi > 5 are conceivable (jBegelman fc Pri nglc 2007). For 
a disc that instead transitions from a standard thin disc to a 
truncated disc with an interior RIAF, a relatively constant 
M is implied at this stage with a decreased radiative effi- 
ciency. Both scenarios of an inner disc remaining extended 
to the ISCO and a truncated disc + RIAF predict different 
M behaviour, which is difficult to distinguish observation- 
ally. 

The phenomenological picture of a variable colour cor- 
rection factor being capable of describing the disc spectral 
evolution observed in BHB outbursts advocates for the de- 
velopment and consideration of models that appeal to an 
evolving vertical disc structure and/or intricate disc-corona- 
jet connection. Recent observational work supports the idea 
of a coupled disc + corona playing a critical role in the 
LHS. A dopting sophistic ated, self-consistent modelling tech- 
niques, iReis et al.l (j2012r i studied the BHB XTE J1650-500 
during its 2001/2002 outburst and found definite correla- 
tions between the reflected X-ray flux and power-law flux 
associated with the corona. These results fit into the inter- 
pretation of a light be nding + collapsing coron a model for 
state transitions (e.g.. iMiniutti fc~ Fabian 20041) and do not 
lend support to the truncated disc picture. An organized 
campaign of contemporaneous and frequent X-ray/radio 
monitoring of Cygnus X-l in the LHS reveal that the re- 
flected ^_X^ay_^ux and the radio flux are positively corre- 
lated (|Miller et alj|2012h. which is consistent with the coro- 
nal ejection model of Beloborodovl (|l999f ). As the promis- 
ing prospect of disc + corona coupling in BHB outbursts 
receives more attention, we hope that the relative measure- 
ments of /coi presented here serve as a useful reference for 
physical models that seek to generate predictive diagnostics 
of spectral hardening in accretion systems. 
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5.2 Physical Range of the Colour Correction 
Factor 

The colour correction factor phenomenologically parame- 
terises our ignorance of spectral alterations arising from 
the highly unknown and complicated combinations of disc 
vertical structure, accretion power dissipation properties, 
magnetisation, and inhomogeneities. Lumping the effects of 
these physical processes into the frequency-independent / co i 
parameter is attractive from the standpoint of limiting the 
parameters in MCD blackbody models; however, whether 
this treatment is an appropriate description of BHB discs in 
all spectral stages remains disputed. Here, we discuss the re- 
sults from theoretical endeavors to constrain the physically 
realisable range of the colour correction factor. 



5.2.1 Radiative Transfer in One-Dimensiona 
Atmospheres 



Disc 



Investigations into the physically reasonable f co \ range and 
the regimes for which the colour-corrected blackbody model 
is justified have focused on one-dimensional atmospheric 
disc models. These atmospheres are assumed to have no az- 
imuthal structure and are constructed at many disc radii 
from self-similar standard thin disc equations. In brief, 
frequency-dependent opacities are computed at each radial 
position, allowing the modification of the locally emitted 
blackbody spectrum to be determined for multiple concen- 
tric annuli. The emergent spectrum, which is modified by 
the vertical structure, is computed for each annulus. The in- 
tegrated disc spectrum is obtained by summing over all an- 
nuli and is fitted with the colour-corrected blackbody model, 
yielding a best-fit f co \. 

Following the technique of self-consistently computing 
the r a diative transfer and ve r tical d isc structure (|Ross et al.l 
ll992l 1. IShimura fe Takaharal (|l995l ) considered the effects of 
electron scattering and free- free emission/absorption on the 
emergent disc spectrum. Regardless of the black hole mass 
and radial position in the disc, the local spectra (i.e., the 
emergent spectrum from an annulus) were adequately fit 
with /coi ~ 1.8-2.0 and / co i ~ 1.7 for M ~ M Ed d and 
M ~ 0.1 MEdd, respectively, suggesting a relatively constant 
/coi over a large range of mass accretion rates expected in 
the HSS. For M < 0.01 M Ec id, which is characteristic of 
the LHS, the local spectrum was not well-described by the 
colour-corrected blackbody model due to ineffective (i.e., un- 
sat urated) Compto ni sation within the disc. 



Merloni et al. (|2000l) extended the study of 



Sh imura fc Takaharal |l995) by considering the effects 
of Doppler blurring, gravitational redshift, and allow- 
ing the dissipated accretion energy to be partitioned 



iretion energ y to be p artitioned 
and disc (|Svensson fc Zdziarskil 



betw een the corona 

119941 1. Assuming all of t he accretion e nergy is dis- 
sipated within the disc, I Merloni et all (|2000l ) found 
/coi ~ 1.8 for M = 0.1 — 0.3 MEdd, fully consistent with 
Sh imura fc Takaharal l|l995l ). Exploring the lower accretion 
rate, M = 0.05 Msdd and allowing a signi ficant fraction 
of the energy to be dissipated in the corona, iMerloni et al.l 
(2000) instead found a variable / co i ranging from ~ 1.9 — 2.7. 
Therefore, if coronal activity is important, / co i cannot be 
considered constant, with the suggested physical range 
being 1.7 < /coi < 3. 



iDavis et all ([2005) presents the most sophisticated 
study to date of the effects of disc vertical structure on 
the emergent spectrum, incorporating bound-free opacities 
and a fully relativistic radiative transfer treatment. For disc 
luminosities, Ldisc = 0.01 — 0.3 LEdd, which spans t he lu - 
minosity regimes of the LHS and HSS, IDavis et al.l (|2005h 
found a weakly evolving / co i ^ 1.4 — 1.6, concluding that 
/coi remains fairly constant over broad ranges in mass ac- 
cretion rate. The results of IDavis et all ( 2005) are at odds 
with those of IMerloni et al.l ( 20001 ); however, a direct com- 
parison is complicated by the different methodologies used 
in each work. 



5.2.2 Magnetisation 

Although the details of magnetic fields (e.g., topology, 
strength, variability) in accreting black hole systems are 
poorly understood, magnetisation is central to many aspects 
of accretion physics. Perhaps most not able is the role of 
the m agnetorotational instability (MRI) (|Balbus fc Hawlevl 
Il99ll ). which is widely accepted as the mechanism for angu- 
lar momentum transport; hence, allowing accretion to per- 
sist. The turbulence generated within the disc by the MRI 
amplifies the toroidal component of magnetic field, perhaps 
reaching a state where magnetic pressu re dominates over the 
combi ned gas and radiation pressures. [Bcgclman fc Pringlel 
(|2007l ) present analytic arguments for the maximal toroidal 
magnetic field strength achievable by the MRI and compute 
the corresponding disc structure, finding that magnetically 
dominated discs are thicker and have harder spectra com- 
pared to the standard Shakura-Sunyaev thin disc. Under the 
conjecture that the toroidal field atta i ns th e limiting value 
set by the MRI, iBegelman fc Pringlel (|2007l ) derive a lower 
limit for the colour correction factor, 



25/108 



/ \ -2/y „ , / \ — 1/lUS / ■ \ b/lUK / \ 

* 6/27 te) (*) (t) 

where p is the mid-plane density, p» is the density corre- 
sponding to the thermalization layer (i.e., effective photo- 
sphere) of the disc, and p*/p < 1. Shearing box simula- 
tions place the disc photosphere a few scale heights above 
the mid-plan e, corresponding to den sity contrasts, p,/p — 
W- 2 - 10~ 4 l|Hirose et al.ll2006l . l2009l l. Therefore, if the disc 
around a stellar mass black hole is magnetically dominated, 
the colour correction factor can far surpass the values pre- 
dicted by thin disc treatments that neglect the effects of a 
magnetic field on the disc structure. 

Conceivably, the details of vertical dissipation and coro- 
nal prope rties could alt e r the emergent disc spectrum con- 
siderably. iBeloborodovl (| 19991 ') envisaged a picture where 
buoyant magnetic loops subject to Parker instability rise 
into the corona and the magnetic energy available within 
the lo ops is subsequen tly dissipated by reconnection. Re- 
cently, Uzdensky] (|2012h analytically computed the vertical 
structure of a gas pressure-dominated disc threaded by a 
poloidal magnetic field where turbulent (i.e., dissipative) 
heating arising from the MRI i s balanced by rad iative cool- 
ing. Based on this disc model. [Uz densky (2012) argue that 
the minimum fraction of available accretion energy that ends 
up being released in the corona is, / m i n ~ Pq^ k ~ 0.4, 
where /3o,max ~ 100 is a reasonable estimate of the gas-to- 
magnetic pressure ratio at the disc mid-plane. Under con- 
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ditions where a substan tial fraction of accre tion power is 
dissipated in the corona, iMerloni et ail ||200Ch demonstrate 
that substantial spectral hardening is expected. 

Numerical simulations show that magnetically sup- 
ported, gas pressure-dominated discs have hard er emergent 
spect r a than their u nmag netised counterparts (iBlaes et al.l 
2006). IBlaes et al.l (|2006h calculated local accretion disc 
spectra from stratified shearing box simulations and found 
the disc spectrum hardened due to magnetic pressure sup- 
port dominating in the upper atmosphere. This magnetic 
pressure acts to vertically extend the disc atmosphere and 
decrease the density at the effective photosphere, leading 
to saturated Comptonisation an d spectral hardenin g. In de- 
riving the local disc spectrum, [Blacs et all (|2006h consid- 
ered two vertical dissipation profiles: (I) a constant dis- 
sipation per unit mass and (2) a broken power-law fit to 
the numerical dissipation profile. Both treatments produced 
identical emergent spectra, implying that the local emer- 
gent disc spectrum is insensitive to the vertical dissipation 
profile and that magnetic pressure support is solely respon- 
sible for the observed spectral hardening. Motivated by the 
result that magnetic pressure support drives the spectral 
hardening and a ssuming the dissip ation profile has no ra- 
dial dependence, IBlaes et all (|2006T ) constructed global disc 
spectra and found only modest hardening, with / co i increas- 
ing from 1.48 to 1.74 upon introducing magnetic pressure 
support. IBlaes et al.l l|2006h caution that the numerical dis- 
sipation occurs at the grid scale and may be sensitive to grid 
resolution. Recent work indicates that numerical dissipation 
in ideal magnetohydrodynamic grid-based simulations does 
behave like physical dissipation, but acts to effectively di- 
minish the resolution (Salvesen et al. 2013, in preparation). 
Adopting dissipation profiles motivated by numerical simu- 
lations is an important step toward merging accretion disc 
theory with simulations. 



5.2.3 Inhomogeneous Discs 

Strong disc inhomogeneities are expected to arise in 
moderate-to-strongly magnetised, radiation pressure - 
dominated regions of accretion discs jBe gclman 2002. [20061') 
due t o photon b ubble instability l|Aronsl 19921 ; iGammid 
1 19981 : iBeeelmanl l200ll ). In this picture, strong density 
contrasts arise on scales smaller than the disc scale height, 
resulting in a porous disc where radiation preferentially 
'leaks' through the underdense regions. Due to the typical 
densities in photon bubble-dominated discs being higher 
than in standard thin discs of the same luminosity, ther- 
malization of radiation is more effective in inhomogeneous 
discs. Consequently, radiation that is thermalized near the 
disc mid-plane and escapes through optically thin chutes 
can result in a colo ur correction fa c tor near unity for 
inhom ogeneous discs (|Begelmanl 120061 ') . iDexter fc Quataertl 
(2012) com puted inhomogeneou s disc spectra based on a 
toy model ^Dexter fc Agoj |201ll) and found that the disc 
spectrum broadens and hardens as the amplitude of the 
local temperature fluctuations increases. Inhomogeneous 
discs can substantially alter the standard disc spectrum 
and further investigation is required to characterise their 
properties. 



5.3 Limitations of Simplistic Accretion Disc 
Models 

While significant strides have been made in exploring the 
imprints of the disc structure on the emergent spectrum in 
the HSS, the applicability of the colour-corrected blackbody 
model to discs in the LHS remains elusive. From a theoreti- 
cal perspective, thin discs with low mass accretion rates will 
be fully gas pressure-dominated with absorption overtaking 
electron scattering as the dominant opacity source, nullify- 
ing the approximation of t he emergent spectrum as a colour- 
corrected blackbody (e.g.. lEbisuzaki et al.lll984T ). Neverthe- 
less, standard MCD models consistently produce exceptional 
fits to LHS accretion disc spectra, giving credence to the no- 
tion that the MCD model successfully describes real astro- 
physical discs in the LHS. Analogously , standard thin discs 
are e xpected to s uffer from thermal ifshakura fc Sunvaevl 
Il976l ) and viscous (|Lightman fc Eardlevlll974l ) instabilities, 
yet the MCD model provides an acceptable description of 
the observed disc spectrum, which should not be the case 
if the disc were plagued by catastrophic instabilities. While 
it may turn out that / co i has a very narrow allowable range 
and/or it is an inappropriate parameterisation to use for low 
disc luminosities, there is currently no compelling evidence 
that this should be so based on fits to observed disc spectra 
alone. 

Direct measurements of / co i should not be made with 
Newtonian disc models when the black hole is known to pos- 
sess a nonzero spin and an inner disc extending to the ISCO. 
Newtonian disc models have the strength of being able to 
describe a standard black hole accretion disc spectrum with 
very few parameters; however, the important physical effects 
of light bending, Doppler blurring, and gravitational red- 
shift inherent to Kerr black holes are neglected. The pseudo- 
Newtonian disc model, diskpn, provides accurate fits to fully 
relativistic disc spectra from a non-spinning black hole, but 
requires uncomfortably large values of / co i when fitting a 
disc spectrum from an extreme Kerr black hole (see JOJ. 
When attempting to fit the hotter disc spectra intrinsic to 
spinning black holes, Newtonian disc models are forced to 
artificially increase the colour correction factor to achieve 
the colour temperature demanded by the data. Fortunately, 
relative measurements of f co \ do not appear to be appre- 
ciably affected in the way that absolute measurements are, 
at least for the low disc inclination considered in this work. 
This implies that the relative measurements of / co i evolution 
during state transitions of GX 339-4 are sound, regardless 
of the intrinsic spin possessed by GX 339-4. The possibility 
that the inner disc recedes radially outward or progresses 
inward between successive observations cannot be ruled out 
and would manifest itself in the relative measurements of 
/coi, rendering them moot. 

The main outcome of this work is that an evolving 
colour correction factor in a standard thin disc with inner 
radius fixed at the ISCO is a viable, albeit not the only, ex- 
planation for the changing disc spectrum of the BHB GX 
339-4 over many state transitions. Given the uncertainty in 
accretion flow geometry as BHBs transition between spec- 
tral states, we emphasise that inner disc radii cannot be re- 
liably computed from spectral fits made with simplified disc 
models. Instead, self-consistent, physically motivated spec- 
tral modelling applied to high-quality CCD X-ray spectra 
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is necessary to distinguish between opposing accretion sce- 
narios. Future insights into the complicated disc-corona-jet 
connections in BHBs will likely progress by concentrated ob- 
serving campaigns and the development of physical models 
that can provide testable diagnostics. 



6 SUMMARY & CONCLUSIONS 

The purpose of this work is to show that physically rea- 
sonable changes in the colour correction factor provide a 
valid, alternative explanation to disc truncation for explain- 
ing the accretion disc spectral evolution observed in black 
hole binary state transitions. To this end, we systematically 
performed spectral fits to RXTE observations of the low- 
mass X-ray binary GX 339-4 over seven state transitions and 
collected the observations that statistically required a disc 
component. The disc detection pipeline measured discs with 
luminosities as low as Ldisc — 10 -3 LEdd and disc fractions 
down to DF ~ 0.05. In a suite of spectral fits to the observa- 
tions deemed to require a disc, we imposed an unchanging 
inner disc radius located at the ISCO for a Schwarzschild 
black hole in order to ensure that a changing colour correc- 
tion factor, rather than a radially receding/advancing inner 
disc, was responsible for any spectral evolution of the disc 
component. We give our conclusions as a bulleted list below. 

• For nearly all state transitions of GX 339-4 investigated 
here, relative changes in / co i in the range 1.4 ^ /coi,rci ^5 5, 
which corresponds to a factor of ~ 3.5 change, can fully 
account for disc evolution in a state transition without 
invoking a truncated disc geometry. 

• Accounting for an / co i evolution in the range quoted above 
over the course of state transitions in GX 339-4, the 
Fdiac oc T G 4 ff relation is recovered for a constant emitting 
area (i.e., non-truncated disc) down to the low luminosi- 
ties characteristic of the HIS and LHS. 

• The validity of the results mentioned above are insen- 
sitive to the Comptonisation model adopted for fitting 
the high-energy continuum. Adopting a continuum model 
that avoids the low-energy flux stealing character of a 
simple power-law reduces the degree of / co i variability 
required to explain the observed disc evolution in state 
transitions. 

• The practice of measuring the inner disc radius from the 
best-fit disc normalisation is unreliable when both R[ n and 
/coi are worked into the disc normalisation. Instead, we 
suggest that disc models where these two parameters are 
decoupled (e.g., diskpn) be used to break the degeneracy. 

• We stress that this work does not demonstrate that the 
inner accretion disc remains at the radius of marginal sta- 
bility during a state transition. Instead, we argue that 
variations of / co i within a theoretically reasonable range 
can account for most, if not all, of the disc spectral evo- 
lution observed over the course of a state transition. 

• Relative measurements, unlike absolute measurements, of 
/coi are not appreciably affected by our ignorance of black 
hole spin. 

• When a physically motivated Comptonisation model is 
used in conjunction with a disc component to fit the 
RXTE data of GX 339-4 in the HIS and LHS, we ob- 
serve a parameter degeneracy between the coronal opti- 
cal depth and the colour correction factor. This calls into 



question the association of ULXs with super-Eddington 
stellar mass black hole binaries based on measurements of 
coronal optical depth without allowing for a variable / co i- 

Simple disc models of the MCD family are attrac- 
tive and useful from the standpoint that they consistently 
achieve superb fits with limited free parameters to accre- 
tion disc spectra, where the intrinsic system parameters are 
often poorly constrained. This work illustrates that allow- 
ing for the possibility of reasonable variations in the colour 
correction factor provides an adequate description of disc 
spectra over a wide range in luminosity states. We advocate 
that /coi may be dynamic and should not be dismissed as a 
mere constant in studies of accretion disc spectra spanning 
a broad range in luminosity, as is often done in practice. 
Numerous factors may contribute to altering the emergent 
disc spectrum including the vertical disc structure, details of 
accretion energy dissipation, coronal dynamics, magnetisa- 
tion, and disc inhomogeneities. Given that the available X- 
ray data cannot isolate these effects, common practice is to 
adopt a blanket term, the colour correction factor, which is 
supposed to account for the overall deviations from a multi- 
colour blackbody disc spectrum. We encourage an increased 
awareness for the implications of presuming a constant / co i. 
Important physics are wrapped up into / co i that should be 
considered before making definitive claims based on the as- 
sumption of an invariant / co i. 

The standard theoretical emergent disk spectrum 
hinges on the assumptions that the accretion energy is dis- 
sipated within the disk (i.e., where the majority of the mass 
resides) and that the density structure is smooth and ho- 
mogeneous. The effects of magnetisation, which surely play 
a crucial role in accretion physics and disc structure, are 
also ignored in the disc models applied to observations. Ra- 
diation pressure-dominated, magnetised discs, as expected 
in the HSS, are subject to clumping instabilities, result- 
ing in a po rous disc struct ure and lowered colour correc- 
tion factor (|Begelmanll2006| y Considerations of the corona 
as a n energy dissipation reservoir (ISvensson fc Zdziarskil 



as a n energy d issipation reservoir (avcnsson &i Z,dziarski 
Il994l ; iMerloni et all |2000| ; [Merloni fc Fabianl 120021) a nd of 
strongly magnetised discs ( Begelman fc Pringld 20071 ') indi- 
cate that /coi may extend well beyond the canonical value of 
1.7 and into a realm consistent with the values required 



here and in similar obse rvational studies (|Dunn et al.ll201ll ; 
iRevnolds fc Miller! [20T1I ') . With the advent of sophisticated 
codes capable of treating magnetohydrodynamics with ra- 
diative transport and general relativity, state-of-the-art nu- 
merical simulations provide a promising avenue for studying 
accretion disc structure with the aforementioned assump- 
tions relaxed. Striving toward the development of sophisti- 
cated models that capture the relevant physics of realistic 
accreting black hole systems, provide predictive diagnostics, 
and are accessible to observers is a crucial step toward disen- 
tangling the details of the accretion flow and distinguishing 
between different accretion scenarios. 
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70 








81(2) 


4.9(2) 


1.26(6) 


0.045(4) 


1.37(8) 


33(8) 


0.739(8) 


70 




16 


-8 


60(2) 


6.9(4) 


1.2(1) 


0.043(8) 


1.2(2) 


46(29) 


0.79(2) 


70 




79 


-8 


65(7) 


6.3(8) 


1.11(8) 


0.038(4) 


1.2(1) 


40(15) 


0.88(2) 


70 


1 


07 


-8 


74(4) 


5.7(5) 


1.05(6) 


0.032(3) 


1.24(9) 


29(10) 


0.914(6) 


70 




98 


-8 


83(2) 


4.7(2) 


0.99(6) 


0.031(2) 


1.20(8) 


35(10) 


0.901(6) 


71 


1 


31 


-8 


82(2) 


4.7(2) 


1.04(6) 


0.033(3) 


1.25(8) 


30(9) 


0.901(6) 


70 




07 


-8 


82(2) 


4.8(2) 


1.05(6) 


0.036(3) 


1.20(8) 


38(12) 


0.901(8) 


70 


1 


19 


-8 


79(3) 


5.0(3) 


1.10(6) 


0.037(3) 


1.28(9) 


29(9) 


0.93(1) 


70 


o 


88 


-8 


79(3) 


5.0(3) 


1.05(5) 


0.036(3) 


1.27(7) 


31(8) 


0.925(9) 


70 


1 


05 


-8 


76(3) 


5.3(3) 


1.01(4) 


0.036(2) 


1.15(7) 


42(13) 


- 


43 


1 


80 


-8 


80(3) 


4.3(3) 


1.10(4) 


0.043(2) 


0.96(3) 


243(36) 


0.96(2) 


70 


1 


16 


-8 


35(2) 


11.5(5) 


1.11(3) 


0.040(2) 


0.867(8) 


1507(81) 


0.99(1) 


70 


1 


76 


-7.742(8) 


26.3(2) 


1.12(3) 


0.041(2) 


0.877(8) 


1431(73) 


0.972(7) 


70 


1 


78 


-7.745(7) 


27.8(1) 


1.05(4) 


0.042(2) 


0.871(8) 


1411(76) 


0.97(6) 


70 


1 


11 


-7.763(9) 


26.4(3) 


1.08(4) 


0.036(2) 


0.890(8) 


1342(70) 


0.99(2) 


73 


1 


1 1 


-7.747(8) 


26.53(6) 


1.12(4) 


0.030(2) 


0.867(6) 


2147(80) 


0.99(2) 


73 


1 


03 


-7.589(5) 


41.5(1) 


0.99(3) 


0.027(1) 


0.866(5) 


2507(84) 


1.02(2) 


72 


1 


48 


-7.523(5) 


39.9(4) 


0.98(4) 


0.024(1) 


0.867(5) 


2591(78) 


0.99(4) 


73 


o 


83 


-7.507(4) 


50.4(5) 


1.04(9) 


0.017(2) 


0.864(5) 


2957(103) 


0.95(8) 


72 


o 


66 


-7.457(6) 


70.0(3) 


- 


- 


0.885(6) 


2491(77) 


1.3(2) 


74 


1 


37 


-7. 


488(6) 


53.6(5) 


0.89(9) 


0.016(2) 


0.803(5) 


3967(154) 


- 


43 


1 


22 


-7.457(7) 


69.4(5) 


1.10(5) 


0.015(1) 


0.885(4) 


2669(72) 


0.87(3) 


72 


1 


13 


-7.457(4) 


71.1(2) 


1.08(5) 


0.028(2) 


0.874(6) 


2017(73) 


- 


43 


1 


67 


-7.601(6) 


44.4(3) 


0.4(2) 


0.0009(3) 


0.557(8) 


5467(443) 


- 


43 


* 




-7.96(1) 


66(1) 


0.76(5) 


0.0014(1) 


0.567(3) 


4401(144) 


1.2(3) 


72 




63 


-8.019(6) 


92.31(9) 


- 


- 


0.50(1) 


6782(1224) 


2.9(6) 


74 


1 


72 


-8 


07(3) 


21(3) 


0.82(4) 


0.0033(2) 


0.640(6) 


1429(87) 


1.3(1) 


72 


1 




-8 


30(1) 


71.7(8) 


0.72(6) 


0.0024(2) 


0.624(8) 


1230(100) 


1.1(2) 


72 


1 


12 


-8 


40(1) 


68.3(9) 


0.96(6) 


0.0040(3) 


0.662(9) 


799(72) 


0.7(2) 


72 


1 


28 


-8 


49(1) 


70(1) 


0.4(1) 


0.0020(4) 


0.58(2) 


1712(334) 


- 


43 


1 


30 


-8 


38(2) 


43.1(4) 


0.81(5) 


0.0032(3) 


0.601(8) 


1492(126) 


- 


43 


* 




-8 


39(1) 


67.4(8) 


0.91(5) 


0.0038(4) 


0.64(1) 


878(97) 


1.3(1) 


71 




73 


-8 


50(1) 


42.7(2) 


0.75(7) 


0.0021(3) 


0.572(9) 


2295(217) 


1.4(3) 


72 




36 


-8 


29(1) 


73.1(5) 


0.88(8) 


0.0021(3) 


0.57(1) 


2581(260) 


2.2(4) 


71 


1 


02 


-8 


23(1) 


57.6(6) 


0.93(4) 


0.0035(2) 


0.647(7) 


1096(75) 


1.0(1) 


71 


1 


^2 


-8 


39(1) 


71(1) 


0.90(8) 


0.0035(5) 


0.58(2) 


2212(377) 


- 


42 


1 


53 


-8 


27(2) 


57.7(6) 


0.78(7) 


0.0035(4) 


0.62(1) 


1270(119) 


1.1(2) 


72 


^ 


58 


-8 


40(1) 


66.0(6) 


0.99(5) 


0.0053(4) 


0.71(1) 


408(49) 


0.9(1) 


72 


1 


19 


-8 


67(2) 


48(1) 


0.96(7) 


0.0049(5) 


0.71(2) 


276(42) 


1.2(1) 


72 


1 


2 r 


-8 


83(2) 


37(1) 


1.22(5) 


0.0058(3) 


0.89(2) 


44(6) 


0.68(3) 


71 


1 


13 


-9 


22(2) 


17.4(8) 


1.0(1) 


0.0034(6) 


0.91(5) 


32(10) 


1.09(8) 


70 


o 


72 


-9 


33(5) 


18(1) 


1.0(2) 


0.0026(7) 


1.06(9) 


11(5) 


0.93(7) 


72 


o 


85 


-9 


51(7) 


12(2) 


0.86(9) 


0.0019(2) 


1.06(7) 


8(3) 


0.92(3) 


72 


1 


18 


-9 


67(3) 


8.6(6) 


1.0(2) 


0.0024(6) 


0.75(8) 


42(34) 


0.86(6) 


72 


o 


88 


-9 


54(8) 


11(2) 


0.9(1) 


0.0020(2) 


0.97(5) 


11(3) 


0.91(3) 


72 


o 


84 


-9 


71(4) 


8.2(7) 


1.0(1) 


0.0021(3) 


0.96(5) 


12(4) 


0.97(4) 


70 


1 


09 


-9 


65(4) 


10.6(8) 


0.8(1) 


0.0013(2) 


0.93(7) 


9(4) 


0.96(3) 


70 


o 


89 


-9 


84(5) 


7.6(7) 


1.0(1) 


0.0014(2) 


1.06(9) 


6(2) 


- 


43 


o 


66 


-9 


80(4) 


7.5(7) 


0.7(1) 


0.0010(1) 


0.96(6) 


7(2) 


- 


43 


o 


91 


-9 


92(4) 


6.1(6) 


1.2(1) 


0.0015(2) 


0.95(6) 


8(3) 


0.92(2) 


70 


o 


99 


-1C 


.04(6) 


4.9(5) 


- 


- 


1.8(2) 


0.5(3) 




45 


1 


09 


-8.0(2) 


6(2) 


0.8(1) 


0.0009(2) 


0.92(8) 


7(4) 


0.94(3) 


72 


o 


87 


-9.96(6) 


5.8(8) 


1.0(2) 


0.0008(3) 


1.6(2) 


0.6(4) 


0.94(7) 


69 


o 


96 


-10.1(2) 


6(2) 


- 


- 


2.3(1) 


0.14(4) 




44 


1 


05 


-10.1(1) 


6(1) 


- 


- 


1.8(2) 


0.2(1) 


- 


45 


o 


73 


-10.4(3) 


5(2) 


- 


- 


2.7(7) 


0.02(3) 


0.79(8) 


73 


1 


19 


-10.6(2) 


3(1) 


0.0(2) 


0.0017(3) 


1.3(2) 


4(3) 


0.91(5) 


72 






-9 


62(3) 


4.2(3) 


- 


- 


1.90(8) 


1.9(3) 


- 


45 


1 


21 


-9 


29(8) 


9(1) 


1.0(1) 


0.0045(6) 


1.00(5) 


20(5) 


0.90(2) 


70 






-9 


36(3) 


7.4(5) 


0.8(1) 


0.0033(4) 


1.02(6) 


16(5) 


0.91(2) 


72 




53 


-9 


45(3) 


5.7(4) 


0.9(1) 


0.0039(7) 


1.1(1) 


11(5) 


0.90(3) 


70 





71 


-9 


43(4) 


6.1(5) 


0.8(1) 


0.0034(5) 


1.21(9) 


9(3) 


0.91(1) 


70 





55 


-9 


40(2) 


6.5(2) 


1.0(2) 


0.004(1) 


1.1(1) 


15(8) 


0.94(3) 


70 





65 


-9 


36(6) 


6.8(8) 


0.9(1) 


0.0047(8) 


0.91(9) 


28(16) 


0.92(3) 


72 


1 


05 


-9 


50(9) 


4.2(7) 


0.8(1) 


0.005(1) 


1.1(1) 


13(7) 


0.93(2) 


70 





89 


-9 


32(4) 


5.7(5) 


0.99(9) 


0.0072(8) 


1.04(6) 


25(7) 


0.98(2) 


70 





91 


-9 


20(3) 


7.1(4) 


1.03(7) 


0.010(1) 


1.12(7) 


28(8) 


0.93(2) 


70 


1 


07 


-9 


02(2) 


9.6(4) 


1.05(6) 


0.012(1) 


1.09(5) 


33(7) 


0.91(2) 


70 


1 


11 


-9 


01(2) 


9.1(4) 


1.08(5) 


0.0128(8) 


1.07(4) 


39(6) 


0.90(1) 


70 


1 


44 


-8 


96(2) 


10.2(4) 


1.18(7) 


0.016(1) 


0.96(4) 


89(17) 


0.96(3) 


70 


1 


18 


-8 


79(3) 


15.0(8) 


1.03(6) 


0.018(1) 


0.78(1) 


694(62) 


0.98(4) 


72 


1 


70 


-8 


27(1) 


36.6(8) 


0.87(4) 


0.0123(6) 


0.770(6) 


1509(73) 


1.19(6) 


72 


1 


27 


-7.948(8) 


52.0(2) 



to 



CO 

a 

C/3 

s 

as 



© 

o 
o 
o 
o 

> 



> 

o 
o 
o 

o 
o 
o 

o 
o 
o 



Trans. 


Observation 


MJD 


Best— fit 


r l 


E h 




r 2 










^linc 


kT in 


^disc 


C 


u 


4 


lo g( F disc) 


DF 


ID 




ID 










Model 




(kcV) 










(eV) 






(kcV) 












crg/s/cm 2 


(%) 


R04 


60705 


01 


70 


01 


53234 


587 


pow X (bb+ga) 


2.42(4) 


- 




- 


0.81(9) 


0.88(4) 


0.0111(6) 


0.726(4) 


2243(81) 


1.1(1) 


72 


1 


76 


-7.880(6) 


62.66(8) 


R04 


90704 


01 


03 


00 


53235 


445 


pow X (bb+ga) 


2.33(2) 


- 




- 


0.63(4) 


1.00(3) 


0.0110(4) 


0.752(4) 


1863(71) 


1.04(5) 


72 


1 


81 


-7.900(6) 


67.7(3) 


R04 


90704 


02 


01 


01 


53237 


547 


pow X (bb+ga) 


2.21(3) 


- 




- 


0.30(2) 


1.01(4) 


0.0072(3) 


0.699(4) 


2629(93) 


1.06(8) 


72 


1 


67 


-7.878(6) 


82.0(4) 


R04 


90704 


02 


01 


03 


53237 


682 


pow X (bb+ga) 


2.31(5) 


- 




- 


0.43(6) 


0.87(5) 


0.0079(5) 


0.676(4) 


3229(126) 


0.94(9) 


72 


1 


49 


-7.847(7) 


77.82(5) 


R04 


60705 


01 


71 


00 


53238 


058 


pow X (bb+ga) 


2.30(6) 


- 




- 


0.39(6) 


0.91(5) 


0.0076(5) 


0.683(4) 


3078(108) 


- 


43 


1 


82 


-7.848(6) 


79.26(2) 


R07 


92052 


07 


6 


01 


54118 


975 


bkn X (bb+ga) 


1.47(1) 


37(2) 


2 


12(6) 


0.52(2) 


0.9(2) 





011(3) 


1.2(1) 


28(12) 


0.93(1) 


70 


1 


32 


-8 


91(3) 


7.3(4) 


R07 


92428 


01 


02 


00 


54122 


304 


bkn X (bb+ga) 


1.49(1) 


35.5(9) 


2 


21(3) 


0.64(2) 


0.8(1) 





013(2) 


1.18(9) 


33(12) 


0.918(7) 


70 




69 


-8 


87(2) 


7.2(3) 


R07 


9 2 -1 2 8 


III 










bkn X (bb+ga) 


1.53(2) 


37(2) 


2 


41(9) 


0.84(4) 


C 


.7(1) 





015(3) 


1.2(1) 


36(18) 


0.93(1) 


70 


® 




-8 


79(3) 


7.3(4) 


R07 






01 


01 


54128 


976 


bkn X (bb+ga) 


1.556(8) 


34.2(9) 


2 


27(4) 


0.97(2) 





95(6) 





023(2) 


0.99(6) 


93(29) 


0.936(7) 


70 




*>ft 


-8 


72(3) 


7.7(6) 


R07 


92035 


01 


01 


00 


54129 


476 


bkn X (bb+ga) 


1.558(9) 


37(2) 


2 


33(7) 


0.99(2) 


1 


02(6) 





024(2) 


1.04(6) 


70(21) 


0.92(1) 


70 


1 


00 


-8 


76(3) 


6.9(4) 


R07 


92035 


01 


01 


03 


54130 


140 


bkn X (bb+ga) 


1.531(8) 


32(1) 


2 


26(4) 


0.96(2) 





97(5) 





029(2) 


0.87(3) 


249(56) 


0.921(9) 


70 


1 


l r 


-8 


52(3) 


11.4(7) 


R07 


92035 


01 


01 


02 


54131 


127 


bkn X (bb+ga) 


1.572(8) 


33.3(9) 


2 


35(3) 


1.14(2) 


1 


03(6) 





027(2) 


1.03(6) 


87(27) 


0.923(7) 


70 


1 


89 


-8 


68(3) 


7.6(5) 


R07 


92035 


01 


01 


04 


54132 


109 


bkn X (bb+ga) 


1.577(8) 


32.3(7) 


2 


35(3) 


1.19(3) 





94(7) 





026(2) 


1.08(7) 


71(24) 


0.937(7) 


70 


i 


43 


-8 


69(3) 


7.4(5) 


R07 


9 2 3 5 


01 


02 


00 


54133 


035 


bkn X (bb+ga) 


1.600(7) 


33.7(7) 


2 


41(3) 


1.30(3) 





90(6) 





027(2) 


1.03(7) 


87(30) 


0.924(6) 


70 


i 


55 


-8 


68(3) 


7.2(5) 


R07 


9 203 5 


01 


02 


01 


54133 


953 


bkn X (bb+ga) 


1.594(8) 


33.3(7) 


2 


42(3) 


1.29(3) 









029(3) 


1.06(7) 


80(27) 


0.916(7) 


70 


i 


41 


-8 


66(3) 


7.6(5) 


R07 


9 2035 


01 


02 


02 


54135 


054 


bkn X (bb+ga) 


1.612(8) 


32.9(8) 


2 


41(3) 


1.36(3) 





94(7) 





027(3) 


1.12(8) 


60(21) 


0.900(6) 


70 


i 


11 


-8 


70(3) 


6.8(4) 


R07 


9 2 3 5 


01 


02 


04 


54137 


018 


bkn X (bb+ga) 


1.640(7) 


32.6(7) 


2 


46(3) 


1.60(3) 


1 


03(5) 






1.05(6) 


84(25) 


0.905(6) 


70 




19 


-8 


66(3) 


6.8(4) 


R07 


92035 


01 


02 


08 


54137 


857 


bkn X (bb+ga) 


1.68(1) 


34(1) 


2 


56(7) 


1.76(6) 





96(7) 





033(3) 


1.11(9) 


69(29) 


0.92(1) 


70 


i 


48 


-8 


65(3) 


8.2(6) 


R07 


9 2 3 5 


01 


02 


07 


54138 


846 


bkn X (bb+ga) 


1.702(6) 


30.4(8) 


2 


49(4) 


1.82(3) 





96(4) 





039(2) 


0.85(3) 


325(80) 


0.935(7) 


70 


i 


68 


-8 


44(3) 


11.0(8) 


R07 


9 2 3 5 


01 


02 


06 


54139 


963 


bkn X (bb+ga) 


1.827(7) 


30.8(7) 


2 


64(4) 


2.39(5) 


1 


06(4) 





040(2) 


1.07(5) 


92(24) 


0.912(7) 


70 


i 


28 


-8 


59(2) 


9.7(5) 


R07 


9 2 3 5 


01 


03 


00 


54140 


225 


bkn X (bb+ga) 


1.871(8) 


30.0(9) 


2 


59(4) 


2.59(6) 





99(4) 





038(2) 


1.02(5) 


130(36) 


0.935(8) 


70 


i 


49 


-8 


52(3) 


11.6(7) 


R07 


92035 


01 


03 


01 


5 4 1 4 1 


075 


bkn X (bb+ga) 


1.979(8) 


31(1) 


2 


74(5) 


3.14(7) 





97(4) 





039(2) 


0.97(4) 


211(49) 


0.937(8) 


70 


i 


38 


-8 


41(3) 


15.4(8) 


R07 


92035 


01 


03 


02 


54142 


056 


bkn X (bb+ga) 


2.088(7) 


31(1) 


2 


82(6) 


3.62(7) 





98(3) 





046(2) 


0.81(2) 


942(129) 


0.944(9) 


70 


i 


75 


-8 


07(2) 


19.2(8) 


R07 


92035 


01 


03 


03 


54143 


039 


bkn X (bb+ga) 


2.26(1) 


26(5) 


2 


.7(1) 


5.1(1) 


1 


00(4) 





046(2) 


0.84(2) 


1114(119) 


0.98(3) 


70 


i 


82 


-7 


92(2) 


20.3(5) 


R07 


92428 


01 


04 


00 


54143 


890 


bkn X (bb+ga) 


2.33(1) 


34(2) 


2.9(1) 


5.4(1) 


1 


02(4) 





044(2) 


0.85(1) 


1391(119) 


0.97(1) 


70 


i 


71 


-7 


81(1) 


23.5(4) 


R07 


92428 


01 


04 


02 


54144 


098 


bkn X (bb+ga) 


2.32(1) 


33(3) 


3.1(2) 


5.5(2) 


1 


01(5) 





044(3) 


0.87(1) 


1201(116) 


0.99(2) 


70 


i 


54 


-7 


84(1) 


22.2(4) 


R07 


92428 


01 


04 


03 


54144 


883 


bknX (bb+ga) 


2.36(2) 


22(8) 


2.8(2) 


5.5(3) 


1 


10(5) 





047(3) 


0.86(1) 


1645(129) 


1.04(5) 


70 


i 


72 


-7 


73 ( 1 ) 


27.0(3) 


R07 


92428 


01 


04 


04 


54145 


974 


bkn X (bb+ga) 


2.41(1) 


41(6) 


3.3(5) 


5.8(2) 


1 


01(5) 





039(2) 


0.86(1) 


1770(125) 


0.94(2) 


70 


i 


51 


-7 


68(1) 


27.3(4) 


R07 


92035 


01 


04 


00 


54147 


031 


bkn X (bb+ga) 


2.51(1) 


29(3) 


2.9(1) 


5.4(2) 


1 


03(4) 





036(2) 


0.876(7) 


2628(123) 


1.00(2) 


70 


i 


48 


-7.484(7) 


38.0(3) 


R07 


92035 


01 


04 


01 


54148 


153 


pow X (bb+ga) 


2.47(3) 








2.1(1) 





84(5) 





019(1) 


0.826(5) 


4196(154) 


1.09(4) 


72 


i 


90 


-7.381(6) 


66.5(5) 


R07 


92085 


01 


01 


00 


54150 


761 


pow X (bb+ga) 


2.36(2) 








1.29(6) 





94(4) 





018(1) 


0.842(4) 


3936(115) 


0.99(3) 


72 


i 


49 


-7.377(4) 


77.3(2) 


R07 


92085 


01 


01 


04 


54151 


816 


pow X (bb+ga) 


2.38(2) 








1.86(9) 


1 


03(4) 





024(1) 


0.838(4) 


3806(123) 


0.97(3) 


72 


i 


74 


-7.400(5) 


69.2(3) 


R07 


92085 


01 


01 


05 


54152 


855 


pow X (bb+ga) 


2.37(2) 








1.04(7) 





91(4) 





016(1) 


0.823(4) 


4454(137) 


1.01(4) 


72 


i 


33 


-7.362(5) 


81.5(2) 



D07 


92704 


03 


()(> 


02 






pow X (bb+ga) 


2 


34(7) 






0.30(5) 


0.71(6) 


0.0048(1) 


0.642(7) 


1939 ( 136) 


1.3(2) 


72 


1 


27 


8.16(1) 


71.2(5) 


D07 


92704 


03 




03 


54223 


293 


pow X (bb+ga) 


2 


40(9) 






0.35(9) 


0.73(8) 


0.0044(6) 


0.62(1) 


2273 (248) 


1.4(3) 


72 


1 


14 


-8.14(2) 


68.3(6) 






!)■■ 






54224 


124 


pow X bb 


3 


06(7) 






1.6(3) 






0.53(2) 


4845 (1518) 


1.9(4) 




| 








D07 


92704 






01 


54225 


253 




2 


34(7) 








75(7) 


0038(4) 


0.633 (9) 


1802 (157) 






43 




64 


8 22(1) 


69 4^7) 


D07 


92704 


03 


08 


00 


54226 


106 


pow X (bb+ga) 


2 


33(4) 






0.34(4) 


0.82(4) 


0.0054(3) 


0.668(6) 


1431(81) 


1.0(1) 


72 


1 


82 


-8.220(9) 


65.4(5) 


D07 


92704 


03 


08 


01 


54227 


624 


pow X (bb+ga) 


2 


08(7) 






0.13(3) 


0.86(6) 


0.0042(4) 


0.633(7) 


1843(137) 


1.1(2) 


72 


1 


85 


-8.21(1) 


82(1) 


D07 


92704 


03 


09 


00 


54228 


399 


powX (bb+ga) 


2 


24(6) 






0.22(4) 


0.76(6) 


0.0038(3) 


0.637(8) 


1575(132) 


1.0(1) 


72 


1 


92 


-8.26(1) 


73(1) 


D07 


92704 


03 


10 


00 


54231 


616 


powX (bb+ga) 


2 


22(6) 






0.24(4) 


0.83(5) 


0.0052(4) 


0.665(8) 


1256(97) 


1.0(1) 


72 


1 


60 


-8.29(1) 


70(1) 


D07 


92704 


03 


10 


11 


54232 


682 


pow X (bb+ga) 


2 


16(4) 






0.22(3) 


0.93(4) 


0.0060(3) 


0.672(7) 


1153(72) 


1 


.1(1) 


72 


1 


94 


-8.30(1) 


69.3(9) 


D07 


92704 


03 


11 


00 


54234 


858 


powX (bb+ga) 


1 


98(3) 






0.26(2) 


0.82(6) 


0.0056(4) 


0.77(2) 


214(32) 





90(7) 


72 


1 


41 


-8.80(2) 


35(2) 


D07 


92704 


03 


11 


01 


54235 


797 


pow X (bb+ga) 


1 


85(3) 






0.20(2) 


0.96(7) 


0.0058(5) 


0.80(2) 


135(23) 





91(6) 


72 


1 


67 


-8.92(3) 


30(2) 


D07 


92704 


04 


01 


00 


54236 


383 


pow X (bb+ga) 


1 


90(3) 






0.25(2) 


0.76(7) 


0.0049(5) 


0.84(4) 


81(21) 





94(5) 


72 


1 


12 


-9.06(4) 


22(2) 


D07 


92704 


04 


01 


01 


54236 


455 


pow X (bb+ga) 


1 


80(3) 






0.19(1) 


0.90(6) 


0.0055(4) 


0.82(2) 


113(18) 





94(5) 


72 


1 


30 


-8.97(2) 


27(1) 


D07 


92704 


04 


01 


02 


54236 


523 


pow X (bb+ga) 


1 


80(3) 






0.20(2) 


0.86(8) 


0.0046(5) 


0.96(4) 


43(10) 





89(5) 


72 


1 


12 


-9.11(3) 


21(1) 


D07 


92704 


03 


12 


00 


54236 


597 


pow X (bb+ga) 


1 


83(5) 






0.21(3) 


0.9(2) 


0.0039(9) 


0.97(7) 


36(13) 





92(7) 


72 


1 


11 


-9.28(4) 


14(1) 


D07 


92704 


04 


01 


03 


54237 


298 


pow X (bb+ga) 


1 


78(4) 






0.18(2) 


1.0(1) 


0.0052(7) 


0.89(4) 


64(15) 






43 





81 


-9.06(3) 


23(2) 


D07 


92704 


04 


01 


04 


54237 


365 


pow X (bb+ga) 


1 


78(3) 






0.18(2) 


0.96(7) 


0.0054(5) 


0.90(4) 


62(14) 






43 


1 


14 


-9.07(3) 


22(2) 


D07 


92704 


04 


01 


05 


54237 


433 


powX (bb+ga) 


1 


74(3) 






0.16(1) 


1.00(7) 


0.0053(5) 


0.94(3) 


51(10) 






43 


1 


21 


-9.07(3) 


23(1) 


D07 


92704 


03 


12 


01 


54237 


501 


pow X (bb+ga) 


1 


78(3) 






0.18(1) 


0.93(7) 


0.0050(5) 


0.89(3) 


62(13) 





89(4) 


72 


1 


67 


-9.08(3) 


22(1) 


D07 


92704 


04 


01 


06 


54237 


704 


pow X (bb+ga) 


1 


80(2) 






0.20(1) 


0.87(9) 


0.0039(4) 


0.96(4) 


36(8) 






43 


1 


65 


-9.18(3) 


18(1) 


D07 


92704 


03 


13 


00 


54238 


750 


pow X (bb+ga) 


1 


66(4) 






0.15(2) 


1.1(1) 


0.0039(7) 


0.95(6) 


31(12) 





95(5) 


72 


1 


19 


-9.42(5) 


11(1) 


D07 


92704 


03 


13 


01 


54239 


743 


bkn X (bb+ga) 


1 


57(3) 


39(5) 


2.8(6) 


0.12(1) 


0.8(2) 


0.0025(6) 


1.10(9) 


13(5) 





93(4) 


70 





98 


-9.38(4) 


13.4(9) 


D07 


92704 


03 


13 


02 


54240 


308 


pow X (bb+ga) 


1 


63(3) 






0.14(1) 


0.6(1) 


0.0024(4) 


0.95(8) 


21(11) 





88(4) 


72 





79 


-9.44(5) 


10(1) 


D07 


92704 


03 


14 


00 


54241 


570 


pow X bb 


1 


49(5) 






0.09(1) 






1.58(7) 


2.4(4) 





91(5) 


74 


1 


13 


-9.51(6) 


10(1) 


D07 


92704 


03 


14 


01 


54242 


884 


pow X (bb+ga) 


1 


52(2) 






0.094(5) 


0.7(1) 


0.0018(3) 


1.05(8) 


10(4) 





82(4) 


72 





83 


-9.58(4) 


8.5(7) 


D07 


92704 


03 


14 


02 


54243 


864 


pow X bb 


1 


48(5) 






0.08(1) 






1.53(8) 


2.0(3) 






45 


1 


07 


-9.64(6) 


8(1) 


D07 


92704 


03 


14 


03 


54244 


979 


pow X (bb+ga) 


1 


52(4) 






0.083(4) 


0.0(2) 


0.0008(1) 


1.2(1) 


5(2) 






43 





89 


-9.72(5) 


7.2(8) 


D07 


92704 


04 


02 


00 


54245 


794 


pow X (bb+ga) 


1 


48(2) 






0.070(3) 


0.8(1) 


0.0015(2) 


0.84(4) 


25(7) 






43 


1 


09 


-9.57(4) 


10.6(9) 


D07 


92704 


04 


02 


02 


54245 


892 


powX (bb+ga) 


1 


51(3) 






0.078(6) 


0.8(3) 


0.0012(4) 


1.0(1) 


11(7) 





89(4) 


72 





95 


-9.80(5) 


6.2(6) 


D07 


92704 


03 


15 


00 


54246 


873 


pow X (bb+ga) 


1 


52(4) 






0.074(8) 


0.0(4) 


0.0005(1) 


1.1(2) 


3(3) 






43 





86 


-9.90(7) 


5.5(7) 


D07 


92704 


03 


16 


00 


54247 


923 


bkn X bb 


1 


49(6) 


29(11) 


1.8(3) 


0.06(1) 






1.4(2) 


1.4(5) 


1 


01(9) 


72 


1 


18 


-9.9(1) 


6(1) 


D07 


92704 


03 


17 


00 


54248 


836 


bkn X (bb+ga) 


1 


46(4) 


29(12) 


1.8(3) 


0.055(6) 


0.9(2) 


0.0011(3) 


1.0(1) 


8(6) 





91(8) 


70 





71 


-9.81(7) 


8(1) 


D07 


92704 


03 


18 


00 


54249 




pow X bb 


1 


48(5) 






0.053(7) 






1.5(1) 


1.1(4) 





91(6) 


74 


1 


23 


-8.0(1) 


5(1) 


D07 


92704 


03 


19 


00 


54250 


731 


pow X (bb+ga) 


1 


45(4) 






0.044(5) 


0.2(2) 


0.0004(1) 


1.2(1) 


2(2) 





72(5) 


72 





83 


-9.99(6) 


6.0(7) 


D07 


92704 


04 


02 


01 


54251 


708 


powX (bb+ga) 


1 


47(3) 






0.042(3) 


1.1(2) 


0.0009(2) 


0.96(9) 


6(3) 





87(4) 


72 





82 


-9.99(7) 


6.6(9) 


D07 


92704 


04 


02 


03 


54251 


780 


pow X (bb+ga) 


1 


50(3) 






0.044(4) 


0.8(2) 


0.0006(2) 


1.1(1) 


3(2) 






43 





69 


-10.09(7) 


5.5(8) 


D07 


92704 


04 


02 


04 


54251 


848 


pow X bb 


1 


50(5) 






0.044(6) 






1.6(1) 


0.5(2) 





99(7) 


74 


1 


18 


-10.2(1) 


5(1) 


D07 


92704 


03 


22 


00 


54254 


925 


pow X bb 


1 


57(6) 






0.043(7) 






1.2(3) 


1(2) 






45 


1 


19 


-10.3(2) 


4(1) 


D07 


92704 


03 


24 


00 


54256 


302 


pow X bb 


1 


51(8) 






0.036(7) 






1.4(2) 


0.6(4) 






45 


1 


05 


-10.3(2) 


5(2) 


D07 


92704 


04 


03 


00 


54256 


612 


pow X bb 


1 


48(5) 






0.031(4) 






1.7(1) 


0.38(8) 





89(6) 


74 


1 


06 


-10.2(1) 


6(1) 


R10 


95409 


01 


12 


00 


55281 


592 


pow X (bb+ga) 


1.53(1) 






0.75(2) 


0.87(7) 


0.019(1) 


0.85(4) 


204(50) 






43 


1 


57 


-8.64(3) 


9.8(7) 


R10 


95409 


01 


12 


04 


55286 


738 


powX (bb+ga) 


1.59(1) 






1.00(3) 


0.80(6) 


0.021(2) 


0.78(4) 


349(104) 






43 


1 


85 


-8.56(4) 


10.1(9) 


R10 


95409 


01 


13 


03 


55288 


380 


powX (bb+ga) 


1.575(9) 






0.97(2) 


0.96(5) 


0.026(2) 


0.80(3) 


367(82) 






43 


1 


98 


-8.51(3) 


11.2(8) 


R10 


95409 


01 


13 


00 


55289 


628 


powX (bb+ga) 


1.587(9) 






1.02(2) 


0.95(5) 


0.026(2) 


0.81(3) 


340(78) 






43 


1 


87 


-8.51(3) 


10.9(8) 


R10 


95409 


01 


13 


05 


55292 


796 


powX (bb+ga) 


1.61(1) 






1.18(4) 


0.93(7) 


0.031(2) 


0.82(4) 


346(100) 






43 


1 


49 


-8.49(4) 


10.6(9) 


R10 


95409 


01 


13 


01 


55293.095 


powX (bb+ga) 


1.625(9) 






1.25(3) 


0.80(6) 


0.025(2) 


0.88(4) 


212(63) 






43 


1 


80 


-8.57(4) 


8.8(7) 


R10 


95409 


01 


13 


06 


55291 


129 


powX (bb+ga) 


1.624(9) 






1.24(3) 


0.97(5) 


0.032(2) 


0.85(3) 


264(63) 






43 


1 


33 


-8.54(3) 


9.4(7) 


R10 


95409 


01 


14 


00 


55295.008 


pow X (bb+ga) 


1.66(1) 






1.32(4) 


0.82(6) 


0.026(2) 


0.80(4) 


352(106) 






43 


1 


54 


-8.51(4) 


10.2(9) 



Ob 

o 

a 
a. 

s 

=1 



to 

05 

a 

s 

to 



s 

to 



to 



Trans 


Observation 


MJD 


Best-fit 


Tl 




r 2 ifpow 




^line 


fcT in 




C 


v 


X u 


l°8(Fdi S c) 


DF 


ID 


ID 




Model 




(koV) 




(eV) 




(keV) 










crg/s/cm 2 


(%) 


RIO 


95409-01-14-02 


55297.883 


pow X (bb+ga) 


1.852(8) 


_ 


2.03(4) 


0.96(4) 


0.037(2) 


0.84(3) 


391(73) 


_ 


43 


1.89 


-8.39(2) 


12.5(7) 


RIO 


95409-01-14-06 


55299.787 


powx (bb+ga) 


2.07(1) 




3.01(8) 


0.95(4) 


0.039(2) 


0.77(2) 


1045(153) 


_ 


43 


1.85 


-8.11(2) 


19.7(9) 


RIO 


95409-01-18-04 


55327.044 


powx Cbb+ga) 


2.43(6) 


_ 


2.0(3) 


0.96(7) 


0.028(2) 


0.768(8) 


3461(198) 


_ 


43 


1.61 


-7.59(1) 


57.3(2) 


Dll 


96409-01-02-02 


55573.469 


pow X (bb+ga) 


2.34(8) 


_ 


0.35(8) 


0.76(6) 


0.0057(6) 


0.61(1) 


2662(271) 


_ 


43 


1.96 


-8.10(2) 


71(1) 


Dll 


96409-01-04-02 


55586.502 


powx (bb+ga) 


2.3(1) 


_ 


0.31(8) 


0.8(1) 


0.0044(6) 


0.62(2) 


1765(285) 


_ 


43 


1.86 


-8.27(2) 


65(1) 


Dll 


96409-01-04-07 


55587.503 


pow X (bb+ga) 


2.2(1) 




0.2(1) 


0.8(2) 


0.0038(9) 


0.60(2) 


1640(422) 


_ 


43 


1.08 


-8.34(3) 


68(3) 


Dll 


96409-01-04-08 


55588.554 


pow X (bb+ga) 


2.3(2) 




0.18(8) 


0.8(1) 


0.0031(6) 


0.57(2) 


2307(590) 




43 


1.35 


-8.28(4) 


75(4) 


Dll 


96409-01-05-04 


55590.441 


pow X (bb+ga) 


2.4(1) 




0.3(1) 


0.6(1) 


0.0031(6) 


0.58(3) 


1787(610) 




43 


1.29 


-8.38(4) 


58(3) 


Dll 


96409-01-05-05 


55592.739 


pow X (bb+ga) 


2-1(1) 




0.15(5) 


0.94(9) 


0.0046(6) 


0.63(2) 


1024(196) 




43 


1-93 


-8.47(2) 


68(3) 


Dll 


96409-01-05-03 


55594.902 


pow X (bb+ga) 


2.11(6) 




0.20(3) 


0.75(7) 


0.0040(4) 


0.66(2) 


571(97) 




43 


1.86 


-8.63(2) 


53(2) 


Dll 


96409-01-06-00 


55597.264 


pow X (bb+ga) 


1.94(6) 




0.17(3) 


1.0(1) 


0.0045(6) 


0.76(3) 


163(40) 




43 


1.40 


-8.94(3) 


35(2) 


Dll 


96409-01-06-01 


55598.677 


pow X (bb+ga) 


1.90(5) 




0.14(2) 


0.79(7) 


0.0041(4) 


0.74(2) 


198(34) 




43 


1.24 


-8.89(2) 


41(2) 


Dll 


96409-01-06-02 


55601.891 


powX (bb+ga) 


1.66(4) 




0.09(1) 


1.05(9) 


0.0041(5) 


0.83(3) 


63(14) 




43 


0.72 


-9.19(3) 


25(2) 


Dll 


96409-01-07-00 


55603.992 


pow X (bb+ga) 


1.67(6) 




0-10(2) 


0.9(2) 


0.0020(6) 


0.8(1) 


28(23) 




43 


0.88 


-9.52(9) 


13(3) 


Dll 


96409-01-07-03 


55604.912 


powX (bb+ga) 


1.63(4) 




0.088(9) 


0.7(2) 


0.0014(4) 


0.76(7) 


45(31) 




43 


0.93 


-9.50(7) 


14(2) 


Tabic 3: 


Best-fit parameters and la standard uncertainty 


in parentheses for 


all observations where 


a disc was 


deemed rcqu 


red. The 


spectral models 


employed for the 


disc and continuum 


were diskbb 



and powerlaw (or bknpower), respectively. The Fe Kq line at 6.4 keV is fit with the gauss model. The entire model is absorbed using phabs with iVjj — 5.7 X 10 21 cm 2 (Miller ct al. 2009b) and an 



energy-independent constant is applied to the HEXTE data to account for offset relative to the PCA. A 0.6% systematic error is added to all PCA energy channels. From left to right, the columns 
give the state transition ID, observation ID, MJD for the midpoint of the PCA observation, the best-fit XSPEC model, power-law index, break energy", high-energy power-law index", power-law 
normalisation, Gaussian line width, Gaussian line normalisation, disc colour temperature as returned by diskbb, disc normalisation, constant applied to the HEXTE data b , number of degrees of 
freedom, reduced x 2 > disc flux over the 0.1-10 keV band, and disc fraction. Dashes indicate that a particular parameter was not applicable in the fit. "Applies for fits using the bknpower model; 
^Applies for fits using the HEXTE data. The full version of this table for all fits used in this work is available as an electronic supplement. 



